DYNAMICS Of SLENDER CYLINDRICAL ELEMENTS CONVEYING 
FLUID AND SUBJECTED TO STEADY AXIAL FLOW 


By 

S. R. BHATE 


nE 

BVN 



DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

AUGUST, 1976 



DYNAMICS OF SLENDER CYLINDRICAL ELEMENTS CONVEYING 
FLUID AND SUBJECTED TO STEADY AXIAL FLOW 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


By 

S. R. BHATE 


to the 

DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUE 

JANUARY, 197" 



LIT. 

CEm^AL L jBRAay 

Acc,. N,. ri 


B 


I'^'^Q'- PO- fbHfi-DyAJ 







ACKWO\€LSDGEMBIIT 


I am extremely thankful to Dr. V. Sundararajan for 
suggesting this topic and the inspiration I received froni 
him during early part of this work. I express deep gra- 
titude to^jards Dr. B. Sahay and Dr. h. Sri Dam for useful 
discussions, encouragement and timel;/ g^aidance. 

I acknowledge with thanks the invaluable cooperation 
and technical assistance received from Messrs. M.I!, Singh, 
S.S. Pathak, R.S. Tripathi, G.K. Tiwari, B.L. Sharma, 

Ram Jiyav/an and C.P. Singh during the design, fabrication 
and assembly of the experimental set-up planned foi’ future 
work. I am thanl^ful to all the staff of Central Workshop 
where most of the experim.ental set-up \j3.s fabricated. 

I must express my appreciation to Messrs S.P. Kapoor, 
M.P. Roberts and W. Singh of the Computer Centre, I.I.T. 
Kanpur for their invaluable help during the computational 
work. 

To Mr, l.D. Verma go my thanks for his accurate 
tjrping. 

Thahis are also due to Mr. Lalta Prasad and Ajodhya 
Prasad of the Department of Jfech. Engg. for neatly dupli- 
cating the stencils. 

I am thankful to Mr. B.L. Arora and S.S. Kushwalia 
for their timely help in tracing the figures. 



Finally, the financial assistance received from 
Bhabha Atomic Research Centre (Grant Mo. DAE/MS/ 7 S- 13 ) 
is "being gratefully acknowledged. 


S.R. Bhate 



TABLE OF COI'ITEHTS 


Certificate 
Acknowle dgement 
List of Figures 
Nomenclature 
Synopsis 


iii 

vii 

xi 

xiir 


Chapter I ; 

IIIERCDTJCTION 

1 

1.1 

Introduction to Flow Induced Vibration 
Problem 

-f 

i 

1.2 

Previous Work on the Lynamics of Pines 
Conveying Fluid 

3 

1.3 

Previous Work on the Dynamics of Rods 
Subjected to Axial Flow- 

14 

1.4 

Present Work, Objectives and Scope 

19 

Chapter II i 

THEORY 

24 

II. 1 

Introduction 

24 

II.2 

Formulation 

24 

II.3 

Dimensionless Parameters 

39 

II. 4 

Response and Stability 

41 

Chapter III : 

THE FLUID INTERACTION PAIAMETERS 

57 

III.1 

Introduction 

57 

III, 2 

Normal Drag Coefficient 

57 

III.3 

Skin Friction Coefficient C^ 

62 

I II .4 

Added Mass Coefficient 

64 





vu 

Chapter IV : 

CO^IHJTATIOifi.L APPEQiiCH 

67 

IV.1 

Characteristic Equations and 
Eigenfunctions 

67 

IV. 2 

Galerkin Matrices 

70 

IV. 3 

S^j-sten Jiatrix and Eigenvalues 

71 

Chapter V : 

EESULTS Ai® DISCUSSION 

7^ 

V.1 

Introduction 

7^ 

V.2 

Isolated Pinned-Pinned Cylinder 

75 

V.3 

Pinned-Pinned Cylinder in Cluster 

79 

V.lf 

Cantilever Cylinder 

81 

V.5 

Pinned-Pinned Pipe 

83 

V.6 

Clamped-Clamped Pipe 

86 

V.7 

Comparison Between Cylinders and Pipes 

00 

V.8 

The liechanism of Instability 

00 

00 

V.8.1 

Buckling Instabilities 

88 

V.8. 2 

Flutter Instabilities 

91 

V.8, 3 

Coupled-Mode Flutter 

95 

Chapter VI ; 

COWCLUSIOITS 

126 

VI. 1 

Summary and Pte suits 

126 

VI. 2 

Future Work 

130 

REFEREITCES 


133 



LIST OP FIGUPES 


vii 


Fig . Ho « 


Caption 


Page 


2.1 Idealized arrangement in pa.rallel flow heat 

exahanger. 53 

2.2 (a) Forces and Moments on a section of c^^'llnder 5^ 

(h) Forces on a section of inner fluid, 55 

2.3 Forces on a section due to mean fluid 

pressure. 56 

5*1 Argand diagram of the complex frequencies, 98 
(jj , of the lowest three rcodes of an 


isolated pinned-pirined cylinder, as functions 
of U 2 ? for p = o-t, €.C^- [ ^ o-/d-i j 
Ti ^1-' ^ 

5.2 Argand diagram of the complex frequencies, 99 

Ct) , of the lowest three modes of an 
isolated pinned-pinned cj^^linder, as functions 
of U 2 J for 

c;^.=.K-Y= IT -T -11^-0 

5.3 Argand diagram of the complex frequencies, 100 
(x) , of the lowest three modes of an 
isolated pinned-pinned c^/linder, as functions 

of Up, for 0*1 ^ € G. = ^ = 1 j 

c:^=-h= y=ti =■ r Jo,-, c-r 


5 , A Ar guild - diagram of the complex frequencies, 101 

do, of the lowest three inodes of an 
isolated pinned-pinned cylinder, as functions 
of U 2 } for Tt •= o • i , b = U = 7 -Tr=: T'= s 
l4 — O * 5~ 

Argand diagram of the complex frequencies, 

00 , of the lowest three modes of an 
isolated pinned-pinned cq^linder, as functions 
of U 2 , for 1 , !cC/^ I , 6 , 

c-<=h=-Y = TT =r^ 0 , 


5.5 


102 



inii 


FIk. IIq. 


Caption 


Page 


5.6 Argand diagram of the complex frequencies, 103 

(j.3 5 of the louest three modes of an 

iso3.ated pinned-pinned cylinder, as functions 
of U25 for p> = 0 * ^ ^ =: I , o<- ^ = TT = 

T = o , u, *- 0*5 

5.7 Argand diagrm of the complex frequencies, 1'C^ 

UJ , of the lowest three modes of an 
isolated pinned-pinned cylinder, as functions 

of U2S for p,- o.45t I > d = 

o<.-k=Y = TT = V^o , 

5.8 Arga.nd diagram of the complex frequencies, I05 

(D , of the i.owest three modes of an 

isolated pinned-pinned cjrlinder, as functions 
of U25 for [i ^ o* 4-8 , €“ 4 o , ^ ~/C-\ , 

c< = K = r= TT = T- = o , -U., ^ Q 

5«9 Argand diagram of the complex frequencies, 106 

LK> , of the lowest fliree modes of an 
isolated pinned-pinned C3^1inder, as functions 
of U2, for 

o( ^ k = r = TP ^ r— o , u., ^ -h 

5.10 Argand diagram of the complex frequencies, IO7 
u.> , of the lowest three modes of an 
isolated pinned-pinned cylinder, as functions 
or "^^2 9 for O * 4" , 4 — 4 d ^ 1 ^ 

C = oc=K= Y-- TT-- r=o , Uj- 


5«11 Argand diagram of the complex frequencies, 

W , of the lov/est three modes of an 
isolated pinned-pinned cyriinder, as functions 


108 


O-L \X^ j 

C 


I or 


j 3 = o- 4-8 , 6 C|. - i j 8 ~ 7 - = ' 
- - ■p - o , i I - o* 5 r .. 1 = 


K=TT - V-o 




y— o 


5.12 Argand diagram of the complex frequencies IO9 
00 , of the lowest three mod^es of an 
isolated pinned-pinned c^rlinder, as functions 
of U2J for p. = = 4-0 5 o - /t - l> 

o(-K - y =Tr-- T - 0 y 'wL,--o*5‘ 


5-13 Argand diagram of the complex frequencies 110 
to , of the lowest three modes of a 
pinned-pinned cylinder in cluster, as j , cr 
functions of Up, for o-t;, 46 = 0-25^ H-f- 

4 , ^ = k u ^ = Y = 13 ^ 



IX 


Fig , Mo , Caption Page 

5 . 1 ^ Argand diagram of the complex frequencies, 111 
CO , of the lowest three modes of a pinned- 
pinned cylinder in cluster, as functions 
of U 25 for ^ - 0 .|;r 6 5 ? 'h - \’5ypC- 

- TT-r' = 0 , u^-O 

5.15 Argand diagram of the complex frequencies, 112 

5 of the lowest three modes of a pinned- 
pinned cylinder in a cluster, as functions 
of Up, for |3 = o* I , 6 Cjr = SFy h.=- 

5 . 16 Argand diagram of the complex frequencies, 113 
cO , of the lowest three modes of a pinned- 
pinned cylinder in a cluster, as functions 

of U^, for 

5.17 Argand diagram of the com.pl ex frequencies, IIV 
(Af , of the lowest three modes of a pinned-*- 
pinned cylinder in cluster, as functions 

of Up, for p~ Q. \ 6. - 1 

0< =: TT ^ "P - 0 , u. i O • £■ 

5*1 8 Argand diagram of the complex frequencies, 115 
VJ , of the lowest three modes of a pinned- 
pinned C 3 rlinder in a cluster, as functions 
of Up? for pi = ~ 4-0 f K_- l » S j & = i 

c><^c.P=^TT=-"i“® ; Ui.2. "s — O* 

5.19 Argand diagram of the complex freq_uencies , 1 16 

Cu> , of the lowest three modes of an 
isolated cantilevered cylinder, as functions 
of Up , f or p> = o « 6 7 € \ , /C = 1 , 

c ^ ^ r = IT = r « d., = o ^ s - Cb= ^ 

5.20 Argand diagram cf the complex frequencies, 

u> , of the lowest three modes of an 
isolated cantilever cylinder, as functions 
of Up, for O'S" , ■c'-4-c> 8 = 7 

IT - 


117 



X 


Fig, Ho . Carjtlon Page 

5 *21 Argand diagram of the complex frequencies, 118 
00 , of the lowest three nodes of an 
isolated cantilever cylinder, as functions 
of u^, for ^ = o* \ 

Y - - TT- r -O , - 

5.22 Argand diagram of the complex frequencies, 119 
00 , of the lowest three modes of an 

isolated cantilever C3'linder, as functions 
of Ups for 6 -o* S' j €r - ^ = 

Y- Hr r = O. 

5.23 Av'gand dicigraii: of the complex frequencies, 120 
60 , of the lox'^est three r;'odes of a pinned- 
pinned pipe, as functions of u^ , for 

e^Vo, Vi - y=Ti-- ‘r = u^=o 

5.2^ Argand diagram of the complex frequencies, 121 
(jj , of the lowest three modes of a pinned- 
pinned pipe, as functions of u-, , for _ 

P Y=C, Ui=. 0.5, 

5.25 Argand diagram of the com*plez frequencies, 122 
(jO , of the lowest three modes of a pinned- 
pinned pipe, as . functions of u^ , for 

o6rt = n -V -h'- Ltj.- O , Ho . 

5.26 Argand diagram of the com.plex fruequencies , 123 

UO 5 of the lowest three modes of a pinned- 
pinned pipe, as functions of u-, , for p=o.S_, 

Ho, 0<= tr p - rr = 6 = K-o ,}6-l,U^=+c*5 

5.27 Agrand diagram of the complex frequencies, 124 

00 , of the lowest three modes of a pinned^:, 
pinned pipe, as functions of u^ , foi^ |3 = 0'b 
TlrP - ^-Ho, -0't> 

5.28 Argand diagram of the complex frequencies 

ix) , of the lowest three modes of a clamped- 
clamped pipe, as functions of u-j, for 
T - U.2.- n - Y = O , <c: = 4-6 


125 



EOllSNCIATUEE 


(Note: The subscripts 1 and 2 refer to the quantities 

related to internal and external flow respectively,) 
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: Cross sectional area of the cj^linder. 

: Flow area, 
s Base drag coefficient. 

: Friction drag coefficients. 

i 

I Equal to (V ) Cj^ L 

s Equal to A: ) and (h /7Z ) C^, respectively. 
; Rotational spring constants. 

: Diameter of cylinder. 

: Hydraulic Diameter = h A^/ wetted perimeter. 

: Flexural rigidity of the cylinder. 

: Inviscid hydrod 3 mamic force on the cylinder, 
per unit length. 

: Normal viscous force on cj’^linder, per unit length 
: Steady state pressiire forces on cylinder, per 
unit length. 

s Geometrj- factor, equal to D/Dj^ 

: Referred to as the angle of incidence, 

I Linear spring constants . 

: Unsupported length of cylinders. 

: Hass of cylinder, per unit length. 

: Virtual mass of fluid, per unit length. 

: Moment on a cross section of the cylinder. 



Axial tension. 

Externally applied uniform tension. 


Time 


Mean flow velocity, cross sectional average. 
)imen£ 


Dimensionless flow velocit 3 ?’, equal to 

-ET- J 

Axial coordinate. 


Lateral displacement of cj/'linder. 
Dimensionless internal damping coefficient 
= 5 i /[6(^ + f,A+M)j 
Dimensionless mass ratio equal to 

Dimensionless uniform tension = T L /El 
Dimensionless gravity parameter 
= (^ m — f-jK +■ f i A ) S' ^ ^ ) 


Angle, the gravitational direction makes with 
X axis , radians 

Equal to 0 or 1 accordingly as the cylinder is 
free or not to slide axially/ at x = L 


Equal to +1 or -1 accordingly as the direction 
of the flow is along positive or negative x axis 
Slenderness ratio, equal to L/D 


Dimensionless displacement, equal to y/L 
Dimensionless spring constants equal to 
k^ L^/ El and kj^ L^/EI, respectivel 3 r 
Dimensionless spring constants equal to 
L^/ El and L^/ El, respectively. 
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: Dimensionless internal damping coeffi :ient 
2 Poisson's ratio 

: Dimensionless coordinate equal to x/L 
•" Dimensionless pressure ps.rameter 
= (p^ - P^) L /Z-I 
i Fluid density 
i Dimensionless density ratio 
: Dimensionless tine equal to 

: Coefficient of virtual mass 
» Circular natural frequency 
J Dimensionless natural frequency, equal to 

^ / E I f'rt ^ 


other quantities are defined in the text 



SYNOPSIS 


A generalized mathematical model has been set up 
to describe small transverse vibrations of slender, straight 
and uniform circular cylinders conveying an incompressible 
fluid at steady velocity simultaneously subjected to a 
steady external axial flow of another incompressible fluid. 
The cylinders may be either ' isolated or be part of a cluster 
of identical cylinders, an arrangement commonly found in 
parallel flow heat exchangers, nuclear reactor fuel rod 
assemblies etc. Gravity, pressurization, surface roughness, 
material damping and directions of the flows has been taken 
into account. Equation of motion derived in the present 
study corrects an error in earlier formulations, related 
to the solid- fluid viscous interaction. The hydrodynamic 
and viscous interaction appears in the form of various co- 
efficients related to the virtual mass, skin friction, 
overall friction for the flov; channel, and normal viscous 
drag, which are functions of Reynolds number, geometrj^, 
compressibility. These have been systematically determined 
from the analytical and experimental results available so 
far . 

The mathematical model has been applied to study 
the dynamic stability of cylinders and pipes with suppor- 
ted ends and clami)ed-free ends. The fluid Interaction, 
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specifically Coriolis forces introduce mode coupling and 
make the problem non-self adjoint. Galerkin Method has 
been used to approximate the continuous system by a dis- 
crete one and dimensionless frequencies of the lowest 
three modes have been obtained with increasing flows. It 
has been shoxm that small flow velocities tend to damp the 
free motions in general but at sufficiently high velocities 
hydroelastic instabilities are possible. Both cylinders 
and pipes with supported ends are subject to buckling 
(divergence) in their first mode. The second mode instabi- 
lity, which sets in at higher flow velocities, is oscilla- 
tory (flutter) for cylinders (in general) and buckling for 
pipes if the mass ratio p is small. At slightly higher 
flow velocities coupled-mode flutter has been obtained for 
supported systems in general, the modes Involved being 
first and second if the mass ratio is small and orrLy the 
fundamental mode if mass ratio is high. For cylinders and 
pipes in cluster, the close spacing severely destabilizes 
the system. Virtual mass of the fluid decides the nature 
of the higher mode instabilities. Cantilever cylinders 
with blunt free end are subject to flutter in all the 
three modes, in addition to buckling in the first mode if 
the direction of flow is towards the fixed end. In 
general the external flow in case of pipes and internal 
flow in case of cylinders is shown to have destabilizing 
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effect on the systems, the effectiveness of the destabiliza- 
tion being dependent on the directions of the flows. Various 
instabilities have been analytically explained in terms of 
fluid-solid interaction. 

Mechanism of energy transfer between the solid and 
fluid has been ana.lyticall3' shown to be responsible for the 
existence of various instabilities. Role of solid fluid 
interaction, both hydrodynamic and viscous, in causing 
various instabilities has been analytically examined using 
energy method. 

Finally a method for obtaining the subcritJcal 
response when the system is subjected to an arbitrary force 
field has been presented. In particular when the excitation 
originates from turbiilence induced random pressiire fluctua- 
tions, the response problem is considered in detail. An 
experimental investigation is in progress to look into the 
type of excitation for a general cluster and to generate the 
data needed to study the response . 
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I . 1 Introduction to Floy Induced Vibratio n Problem 

Hie term "flov/ induced vibraticn-’ includes a wide 
range of vibration problems associated witn flexible solid 
structures activated by the energy/ derived from the fluid 
flowing around or through the structural rembers . High 
velocity coolant flowing through a reactor core is a source ■ 
of energy that can induce and sustain vibration in reactor 
core components. Various components of nuclear reactor, 
such as fuel pins, control rods and heat exchanger tubes, 
are long and slender v/iiich may be set to vibration by the 
forces arising from the fluid flow. Structural vibrations 
often cause damage to such components through wear, fatigue 
and fretting. In fact, there have been numerous flow 
induced vibration problems, not3.bly in EER-II [ 75 3 j 
GETRi; 76-77 3 , Big Rock Point C 78-81 3 , Yankee 182-8^3 , 
SEERT - III [ 77-,85 1 , and FiSPE [86, 873 reactors . 

These problems included i) sustained oscillations, 
which resulted in excessive wear, fretting a.nd fatigue, 
ii) incorrect flow distribution, iii) broken specimen 
holders, iv) excessive vibration in heat exchangers and 
v) fluctuations in the pov/er level of the reactor. The 
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examples of flow induced vibrations in other fields 
are numerous. Critical loads are Imposed on struc- 
tures by gusts and large turbulent motions in the 
atmosphere where the length scale of the turbulence is 
comparable to the size of the structures. The high 
frequency pressure fluctuations of small scale turbu- 
lence can cause a fatigue fail’ure in a surface exposed 
to flow by providing excitation near a resonant fre- 
quency of the structure. A somewhat similar mechanism 
plays a role in the generation of water waves | resonant 
waves are excited b3^ the pressure pattern convected 
over the surface in turbulent wind. Flow induced 
fluctuations can excite aircraft fuselage pannels 
causing noise transmission. 

Vibration of elastic structures due to flue-' 
tuating lift forces set up as a result of vortex 
shedding process in the cross flow of a fluid have been 
observed since a long time due to common occurence of 
this phenomenon in nature . The vibrations due to 
parallel flow were observed much later, only after 
fluid carrying conduits, especially flexible ones, came 
into use,. Indeed, the peculiar spontaneous motion of 
the free end of a flexible pipe by a sufficiently high 
flow rate was recognized as a self excited oscillation 
by Ifercel Brillouin in 1885, but his work on the sub- 
ject remained unpublished. 
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Analysis of the vihratlon prohlems may he diirided 
into the considerations of the stability and response 
problems. In stability analysis investigation is made 
to determine whether conditions exist that can canse the 
structural response to diverge and, if so, to find these 
conditions (stability limits) in terms of the system para- 
meters. In response anal^^sis the problem is to predict 
the actual motion of the structure . The solution of the 
stability problem often supplies sufficient information 
to consider the vibration problem as solved. If the 
actual motion of the structure is desired, or if one 
wishes to know wliat happens to the structure under unstable 
conditions, the response problem must be studied in detail. 
This requires the knowledge of the damping mechanism and 
flow originated excitation wluch may be due to unsteady 
pressure field or piilsations in the flow rate. The res- 
ponse studies of flow induced vibration often necessitate 
solving a nonlinear differential equation coupled with 
the flow dynamics, the analysis might then be termed as 
"fluid-elastic analysis". 

1.2 Previous Work on the Dynamics of 
Pices Conveying fluid 

As stated earlier, the para3.1el flow induced 
motion was first observed in 1885 by Brillouin. 
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Bourrieres Ll3«one of his students, wa,s first to undertake 
the study of the d 3 ’'nainics of fieri 131.9 pipes conveying 
fluid. He examined the oscillatory ins tahili ties of 
cantilever pipes conve 3 /ing fluid goth theoreticaJlly and 
experimentally. He derived the correct equation of motion 
and, although unable to obtain anai.^rticallj’ the critical 
velocity for the onset of oscillations, he determined 
most of the salient characteristics of the phenomenon. 

Interest in the subject uas reactivated in 
19^05 in connection with the study of vibrations of 
Trans-Arabian pipe line, by Ashley and Haviland 12 3. 

They studied the flow induced damping in simply suppor- 
ted pipe and found it to increase rapidly with the flow 
rate. They further studied the >/ave type solutions and 
found that the effect of the flow is to expedite the '&/aves 
travelling in the flov; direction and retard the waves 
travelling in opposite direction. The same problem was 
studied independently by Eousner\.33 using a different 
approach. He found that for sufficiently high flow 
velocities the pipe may buclLLe, essentially like a column 
subjected to axial loading. A subsequent elegant and 
more general study by Niordson 1.63 led to the same equa- 
tion of motion and to essentially the satie conclusions 
regarding stability of pipe with simply supported ends. 
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Long ^ s'} ms the first, after Bourrieres, 
to he concerned with cantilever pipes conveying fluid, 
although he considered other boundary conditions. 

His method of solution, using a power series for the 
mode shape, ms applicable to relatively small flow 
velocities, considerably below the threshold of 
oscillatory^ instability, of the existence of v/hich he 
seemed to be unaware. Nevertheless he did include the 
centrifugal force term in the equation of motion and 
confirmed experimentally that in contrast with those 
of siraply supported pipes, forced motions, of canti- 
lever pipes are damped by internal flow in the range 
of flow velocities . The conclusion of Ashley and 
Haviland \_2'} that the flow of fluid damps the motion 
of the simply supported pipes was due to improper 
selection of lesser number of terms of the power series . 
The damping associated x\dth cantilever pipe originates 
from Coriolis force- and asymmetric modes about the 
centre of the span. 

Hands Iman \_73 presented an analytical method 
in which the character of the eigenvalues of the problem 
is determined from the structure of the differential 
equation of motion without determining specific solu- 
tions. Using a perturbation series in terms of flow 
velocity for the mode shape and eigenvalues he derived 
the same conclusion as Long \_8 3 in case of pipes with 
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supported ends that no damping takes place x^hen pertur- 
hation terms up to second and higher order are retained. 
For fixed-free ends Coriolis force promotes damping for 
small flow velocities. For near critical velocities 
the frequency of vibration is real for subcritical velO' 
cities and purely imaginary for velocities above criti- 
cal; for fixed-free ends, however, there is no critical 
velocity for which frequency of vibration vanishes. 

In all the above studies, excepting 
Bourrieres, the only form of instability discovered >ra.s 
buckling. It was not until 1963 that Gregory and 
Paidoussis ^26, 27^ showed theoretically and experimen- 
tally that at sufficiently high flow velocities, canti- 
lever pipes are subject to oscillato3ry instabilities 
(flutter) rather than buckling (divergence). The 
effect of floiiT for small values of flow velocity is to 
damp .all the rfndes. For higher values of flovr velo- 
city, however, the second mode becomes less dam.ped and 
subsequently the pipe becomes unstable. The buckling 
type of instability does' not occur. Internal damping 
of the material >ra.s found to destabilize the pipe. The 
existence of oscillatory instabilities was fully anti- 
cipated b 3 '- Benjamin £l1 , 12^ in a study concerning the 
dynamics of articulated pipes which is a discrete- 
representation of the continuous flexible system. 
Benjamin was the first to perceive the dynamical 
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problem independent of the fluid friction and he analy- 
tically predicted the existence of oscillatory^ 
instability of cantilever pipes conveying fluid. Both 
these effects later were confirmed by Gregory and 
Paidoussis '[^26, 27j . The dymanic behaviour of arti- 
culated and continuously flexible pipes conveying fluid 
is not strictly analogous. Benjamin [l1, 12"^ found 
that the buckling instability is possible in the case 
of a vertical cantilever articulated system, where 
gravity, is operative, if the fluid has sufficiently?' 
high density; on the other hand, Paidoussis (^37 1 found 
that vertical continuous flexible pipes are never sub- 
ject to buckling. This paradox was clarified by 
Paidoussis and Deksnis {^38 J . The stability of canti- 
lever pipes was further discussed by Ifemat-IJasser , 

Prasad and Hermann £25 3with emphasis on the effect on 
stability of velocity dependent forces, such as dissi- 
pative and Coriolis forces; they showed that such 
forces may destabilize the sy^stem, an effect also repor- 
ted previously by Gregory and Paidoussis [_26, 27l« 

Subsequent work by Hermann t.283 and Hermann 
and Memat-Nasser [333 stressed the connection between 
the problem of instability of a cantilever conveying 
fluid and the more general problem of instability of a 
cantilever subjected to a ’’follower” type of force at 
the free end, i.e., a force retaining the same angular 



8 

disposition relative to the free end during small 
motions of the cantilever. 

A recent study hy Stein and Torbiner L39 1 
mainly concerned with infinitely long jDipes conveying 
fluid, brought a correction to the equation of motion 
as derived in some of the earlier papers. This correc- 
tion arises from the effect of internal pressure and 
may become significant for sufficientl 3 r high pressure. 
Nagulesxi/aran and Williams [353 studied this aspect of 
the problem both theoretically and experimentally. In 
case of pipes supported at both ends the effect of 
internal pressure was found similar to that of flow 
velocity; the pipe ma,y buclile even at very small flow 
velocity by the action of sufficiently high internal 
pressure . 

Thurman and Mote [363 presented a nonlinear 
analysis for a pipe with simply supported ends conve 3 '-- 
ing fluid, using a perturbation technique. Theq^ found 
that in determining the natural frequencies of the 
system, the importance of the nonlinear terms increases 
with the flow velocity, so that the range of applica- 
bility of the linear theory becomes more restricted 
as the flovf velocity increases. 

More recently Chen studied the stability 

of a pipe conveying fluid, with upstream end clamped 
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and downstream end constrained by a linear spring 5 
so that the boundary conditions are intermediate bet- 
ween clamped-free and clamped-pinned ; accordingly both 
buckling and oscillatory instabilities are possible, 
in general, depending on the spring constant. Both the 
types of instabilities with multiple stable and unsta- 
ble ranges of flox-/ velocity may occur for a range of 
spring constant values . For the value of dimension- 
less spring constant below a certain critical value 
(3^*814-5) the pipe does not loose the stability by 
buckling. 

' Paidoussis and Denise C4-3, 4-63 studied the 

dynamics of very thin elastic pipes conveying fluid, 
by utilizing thin shell theory to describe the motions 
of the pipe and potential flow theory to obtain the 
fluid forces. They analyzed both cantilever pipes and 
pipes with clamped ends. They found that in addition 
to instabilities in the beam modes of the system 
(corresponding to those found previously by beam theory 
for thicker pipes), instabilities in the shell modes 
are also possible, as verified 63 / their experiments. 

Of particular interest was the finding that thin pipes 
with clamped ends are not only subject to buckling but 
also to coupled-mode flutter. For supported pipe as 
well as cantilever pipe the system loses stability by 
flutter in circumferential modes, in addition to usual 
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beam modes. Similar theoretic?! results were obtained 
later by Weaver and Unny {.4-91 in the case of simply 
supported shells. The pipe was found to buckle in the 
beam mode at first but, at higher flow rate, instability 
by coupled mode flutter was predicted. The critical 
velocity for the coupled mode flutter depends on the 
circumferential mode number and the mode number for 
minimum critical velocity considerablj?" depends on the 

i> 

length and thickness ratios. The critical flow velo- 
cities for short, thin shells are associated with a 
large number of circumferential modes; further, the 
longer and/or thicker the shell lower the critical 
circumferential mode number until the instability mode 
is the same as in a simple beam. Weaver and Unny also 
studied the divergence boundaries and critical mode 
shape, both theoretically and experimentally. 

In all the studies discussed above the flow 
velocity was taken to be steady. Fbcently Chen £4-1 3 
examined the stability of simpty supported pipes con- 
veying fluid with a flow velocity, U, which has a 
time dependent harmonic component superposed on the 
steady velocity, such that U = (1 + ^ucos CO t). 

His equation of motion incorporated the initial cur- 
vature, damping and axial tension in addition to 
unsteady flow. Using Galerkin method the governing 
equation was reduced to a system of coupled la thieu- 
Hill type equations with multiharmonic coefficients. 
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He investigated the onset of parametric resonance using 
Hsu’s C 663 and Bolotin’s C 65] method and obtaiiaed the 
boundaries of the first two instability regions ,* more- 
over, he found thci.t the combination of resonances is 
also possible. He studied the significance and. effects 
of Coriolis force on free vibration, forced vibration 
and parametric resonance. The influence of Coriolis 
force on free vibration and parametric resonance was 
found small for a flow rate much less relative to the 

critical velocity. Initial curvature, which introduces 
in 

nonhomogeneity/ the differential equation of motion, 
may induce large deflections if the flow velocity is 
large. This particular aspect of the problem was stu- 
died by him in detail ^^2, Mf 3 for a uniformiy curved 
tube conveying fluid.. He included the effect of fluid 
pressure in the equation of motion. Of particular 
interest was the finding that the curved, tube with 
supported end loses stability at first by buckling in 
the first mode, then at higher flow velocity again by 
buckling in the second mode and so on. The flutter 
type of instability was not found to occur in case of 
tubes with supported ends. The effect of the fluid 
pressure on natural frequencies of the tube was found 
quite similar to that of flow velocity. This was already 
found for the straight pipes by Stein and Torbiner [39 3 
and Naguleswaran and Williams [353 • Coriolis force 



and internal friction did not affect tlie stability of 
curved tube with supported ends, as found in the case 
of straight tube Cl1 , 12 3 . Further ^-ork of GhenC^tJ 
dealt with out of plane motion of a uniforml^r curved 
pipe. He derived essentiall;’- same conclusions as those 
for the in-plane motion; thus the pipe experiences 
bucliling type of instability for conservative cases 
(supported ends) and flutter type of instability for 
nonconservative (cantilever pipe) cases as the flow 
velocity is gradually increased. No coupling between 
in-plane and out of plane motion was found and as such 
these can be studied independently. 

Chen's work on parametric instabilities of 
simply . supported pipe was followed by the work of Bohn 
and Hermann C 4-7 3 wherein they examined the transverse 
motion of vertical articulated pipes when the flow rate 
harmonically varies about a nonzero mean. Both parame- 
tric and combination resonances were found to be 
possible. The amplitude of the pulsations in the flow 
rate required to cause the insta.bility was found to 
decrease as the mean flow velocity approached the one 
required to cause the onset of instability either by 
flutter or divergence when the flow is stead^r.. For the 
values of mean flow velocity greater than the one which 
causes t];^' onset of instability by flutter or divergence 
when the flow is steady, there exists a range of 



pulsation frequency and amplitude which can stabilize 
the system, where in the absence of pulsations the 
system would be unstable* Some of the regions on 
Strutt diagrams (curves of pulsation frequency/ versus 
amplitude showing stable and unstable zones) which 
constitute stable regions when mean flox/ velocity is 
less than critical become the regions of combination 
resonance when mean flow velocity exceeds the critical 
value . 


In a recent paper, Paidoussis and IssidC553 
considered the dynamic stability of pipes conveying 
fluid, where the floi-?' velocitj'" is either entirely/ 
steady or with a small harmonic component superimposed 
on it. He rederived the pertinent equation of motion 
correcting an error in Chen’s formulation. This correc- 
tion arises from the streamwise local acceleration 
( -:=:-r ) of the fluid in unsteady flow and had not been 
considered by Chen. Of particular interest was the 
finding that conservative systems were subject not onl 3 r 
to buckling at sufficiently high velocities but also 
to flutter at higher flow velocities. The existence of 
these oscillatory instabilities in case of pinned-pinned 


and clamped-clamped pipes ms explained as the conse- 
quence of Coriolis forces, the system being srroscopic 
conservative., Accordinglj’’ the method of equilibrium 
falls to give correct higher critical velocities. 
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After the onset of buckling the Coriolis forces 
stabilize the systen prior to the onset of co’apled 
rode flutter. He also consider'^d the case of canti- 
lever pipe 1-ath viscoelastic and hysteretic damping. 

The viscoelastic damping ves found to destabilize the 
pipe whereas the effect of h3^steretic dissipation was 
not as severe. For the case of pulsatile flow velocity*, 
primary and secondaipo principal instabilities were 
investigated, using Bolotin's method and stabilitj^ maps 
were presented for pipes with clariped or pinned ends 
as well as cantilever pipes. The extent of the insta- 
bility region was found to increase with flo’/ velocity?' 
for pipes with supported ends, while a more complex 
behaviour was found in the case of cantilever pipes. 

In all cases, dissipation reduced (or entirely elimi- 
nated) the extent of parametric instability zones. 

1,3 Previo us V/ork on the Dy n amics of R ods 
Subjected' to Axi al Flow ; 

That the vibration caused by the energ^g deri- 
ved from the parallel floi’ of a fluid along a cj^lindri- 
cal bod3^ could be significant was not evidenced until 
1958 when Eurgreen et. al. C 9 1 de^'onstrated the pheno- 
menon of rod vibration due to para.llel flow, while 
studying the heat transfer of water flowing across a 
bundle of rods simulating a nuclear reactor core. This 
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was perhaps the earliest study notivated ty the general 
concern that parallel flow induced vihratlon of cj-^lin- 
drical fuel elements may cause a niechanicc'l design 
problem in reactor. Oscillations were observed to occur 
at all flow velocities. The fact that the frequency of 
oscillations remained relativels^ constant led them to 
believe that the vibration was of self-excited type. 

They developed an empirical expression for maximum 
amplitude of vibration by dimensional analysis of 
differential equation of motion in which excitation was 
assumed proportional to the dynamic pressure ( ^ f ) 
of the fluid. This was follovred by several experimental 
studies (*18, 19 , 20, 22, 23, 2k, 32, 3^1 wherein the 
principal aim was to study the character of the vibra- 
tions and obtain an empirical expression for maximum 
amplitude of vibration which can be used for design 
purposes. Tho first analytical and experimental investi- 
gation of the dynamics of flexible slender C 3 'linders 
subjected to axial flow wa.s attempted in 1966, about eight 
years after the work of Eurgresn ot . al . £9l bj^ 

Paidoussis |29, 30^ • The analj’^tical model devised for 
this purpose vras based on Bornoulli-Sulor beam theory 
to describe the flexural motions of the cylinder, 
slender body theory for the coupled invi sc id hydro- 
dynamic forces, and fairlj' simple linearized relation- 
ships for the corresponding viscous force s. The 



cylinders were either supported at totii ends or were like 
a cantilever with a suitably streamlined, free downstream 
end. The theoretice.lly determined main characteristics of 
the dynar'iical behaviour of the system were found to be 
the following; i) small flov? velocities have a damping 
effect on free motion of the s3''stem, ii) for sufficiently 
high flow velocities the S3''stem is subject to h3rdroelastic 
instabilities; both divergence and flutter are generalli^ 
possible, the former representing a static loss of 
stability and the latter being an oscillatory instability. 
Paidoussis showed that the C3'"linder supported at both 
ends represents a conservative g^/roscopic system. 
Accordingly, damping effect at low flow velocities occurs 
entirelj^ due to dissipative forces. Divergence was 
dominant in case of cylinders with supported ends though 
flutter also occurs at higher flox/ velocities in higher 
flexural modes. Those oscillatory- instabilities arc 
speeifically. caused by lateral friction forces and in the 
absence of the drag forces only buckling is possible. 
Externally applied tension >ra.s found to stabilize the 
system and compression destabilize it. A cantilever 
cylinder in axial flow is inherently a non-conservative 
system, so that work done "by the inviscid hydrodyrnamic 
forces in the course of free oscillations of the systosi 
need not be zero. Indeed, at small... floxf velocities this 
work is negative resulting in an added damping effect 
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over and above that due to dissipation. At higher flcv 
velocities the worl; done by the in^’^iscid hydrodjrnanic 
forces ma}^ become positive, so that the energy then may 
be gained by the cylinder from the fluid | indeed this is 
the dominant mechanism of energy transfer underlying the 
onset of flutter of cantilever cylinders. 

Recently Paidoussis [ } extended the theory?- 

and applied it to the c^.'linders in confined flows. He 
took into account the gravity and pressurization effects 
which were not considered in his previous work on C3^1in- 
ders I in addition he corrected the relationship for 
normal friction drag, adding a term to account for the 
flow Induced damping at zero flow velocity. Ihe cylin- 
ders could be either Isolated or be part of a cluster 
of identical cylinders. Smsll flow velocities were found 
to damp free vibration but at sufficiently high flow 
velocities possibilities of hjrdroelastic instabilities 
were found. Typicall3^, simply supported cylinders were 
subject to buckling in their first mode, si nd then at 
higher flow velocities in their second inode; at slightly 
higher flow velocities coupled-mode flutter i/as shox<m to 
occur. Critical velocities I'or the cylinder in a closely 
spaced cluster were much lower than those for an isolated 
cylinder. Ca.ntilever cj^-linders, on the other hand, were 
subject to buclcling in their first mode and to flutter 
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in higher modes . The difference in the djmamic behariour 
of the cylinder as compared to the one in iiis previous 
stud;;- ^ 29 ? 30 J was attributed to the correction in the 
representation of fluid frictional forces incorporated in 
his later study. This demonstrated the definite impor- 
tance of fluid frictional forces in the dynairic behaviour 
of cylinders. Coupled-mode fliitter -was shown to occur for 
confined cylinders also; for cantilever cylinders however 
only ordinary flutter insta,bilit;r \/as found to occur, 
i ore over 3 for the case of cantilever cylinder with stream- 
lined free end the loss of stability in the second mode 
preceeded by regaining of stability in the first, and in 
analogous manner, the regaining of stability in the second 
mode occured at a slightly higher value of floi/ velocity 
than was necessary for the onset of -third-mode flutter. 
Apart from the study of dynamics he presented a raethod by 
-which response of the cylinder subjected to an arbitrary 
force field may be obtained. 

Further work of Paidcussis \_6'l 3 deals -'uith the 
stability of slender cylinders subicctecl to axial flow^ 
the velocity being harmonically perturbed in time , 
Cylinders with pinned-pimied ends and clanped-free ends 
were considereci . Frincipa.l primary and secondar;/ insta- 
bilities xrere found to be possible for certain ranges 
of frequenoies and amplitudes of the perturbations. 



For small flow velocities flow induced damping opposes 
the onset of parametric instabilities. The parameters 
which deststbilize the system when flow velocity is steady 
render parametric instability regions wider. In case of 
simply supported c^/lindcr parametric instabilities occ’or 
over specific ranges of flow \’-elocitics . In the case of 
cantilever cjrlinders the instabilities were found to be 
associated v/ith only some of the modes of the system and 
resulting regions in strutt dia.granis were quite complex^ 
moreover, in certain cases the secondary instability 
regions became more important than the primary/ ones. In 
general, however, the extent of the instability regions 
in strutt diagrams became wider as the system approached 
the point of instability in steady flow. 

1.4 Present Work. Objectives a.nd Scope ; 

From the review of previous work on the dynamics 
of pipes conveying fluid and cylindrical rods siit j acted 
to axial flow, it appears that a remarkable development 
has occured since Trans-Arabian pipeline i-/as obs8r‘';'ed to 
vibrate (presumably as a lesult of internal flow) and 
cylindrical rods were observed to oscillate br/ Burgrten 
et. al. However the impetus for these studies, espe- 
cially those related to pipes,, did not conic from desire 
to solve a pra.ctical problem, but rather from one or 



more of the following considerations i i) the physical 
problem was inherently intriguing; ii) the protler of fere 
scope for interesting mathematical manipulations; 
iii) the pipe problem represents one of the few cases 
khere a (non-consGrvative ) follov/er loai is pliysicallrp 
realizable, hence, it offers a rare opportunj.ty for com- 
bined theoretical-experimental studies in the general 
area of dynamics of elastic systems subjected to non- 
conservative forces, and. expla.ins, to some extent, the 
large volume of research acti^tity on this topic. It 
further appears that the development in the subject of 
dynamics of pipe transporting fluid and in the subject of 
slender rods subjected to fluid flow occured more or less 
independently. These two problems have a unique feature 
in common viz. the energy transfer between the flowing 
fluid and body by inviscid hydrodynamic forces. The 
centrifugal force set-up by the fliiid flowing through 
a curved tube or along a curved rod, a.cts much like a 
compressive force at the end. The problem of rod differs 
from that of pipe oiilj?" due to added frictional forces 
which can also be considered for a pipe- submerged in a 
fluid, o.g. undervreiter pipeline. As regards vibrations 
due to pressure; fluctuations in the flow either due to 
turbulence in the flow or due to its pul!. satile nature , 
both pipes and rods vibrate or loose stability in a 
more or less similar manner. The first aim of the 



present work is to s^mthesize these tv;o problems into a 
single general parallel flow problem, so that in appli- 
cations where a rod subjected to parallel flow also acts 
as a tube conrej/'ing fluid, the effect of both the flows 
on the djmamics of the rod can be assessed. 

The present ’.rork was initiated with an aim to 
obtain an estimate of the response of structures compri- 
sing slender cylindrical tubes or rods subjected to 
fluid flov/ with prospective application to parallel 
flow heat exchangers, liquid fuel rocket piping, nuclear 
reactor coolant channels, fuel rod bundles, underwater 
pipelines etc. Accordingly, the equation of motion 
set up taking into account every possible realistic 
effect. Model for viscous interaction between the tube 
and external fluid was chosen and various fluid parameters 
v/ere systematically determined to obtain the correct 
estimate of damping. Hox'/ever, since the mechanical struc- 
ture was regarded as the sj'stem excited by pressure fluc- 
tuations in the flow (internal and external), either due 
to turbulence or pulsatile nature of the flow, the quanti- 
tative knov/ledge of excitation was a must before response 
could be obtained. Except for the case oi a flat pla,te 
^island circula.r pipe 16, 17lno reliable measure- 

ments on pressure fluctuations due to turbulence are 
available so far. Analytically these can be obtained 



Toy nianerical siniulation of the flow field oiily for sone 
simple geometries such as annular reometry ['62 3tiut the 
method is both difficult and conputationally uneccnoEiical 
due to large computer memory core and time requirements. 
Therefore an experiraerxtal investigation has teen planned. 
In the meanwhile with the available information at hand 
a major portion of the solution of the problem could be 
s ought . 

For any dynamical system, for a given excitation 
the response depends on the natural frequencies and asso- 
ciated damping. The vibration problem in most cases, is 
considered solved if tho naturail freouencies and damping 
is known. In addition, if the conditions at which the 
system tri-ll loose the stability are known in advoiiice, thee 
can be avoided at the design stage itself. Accordingly, 
an investigation of the effect of internal and external 
flow on the natural frequencies and damping taking into 
account the gravity, pressure, material damping and 
confinement in case of. clusters was planned. Since fluid 
friction forces play a,n important role in the dynasiics of 
rod subjected to axial flox/, consideration of the surface 
roughness and viscosity of the fluid is important ^'hich 
has not received due attention so far, presumably because 
in all the previous formulations the fluid friction 
appears as skin friction coefficient and normal drag 



coefficient which have to he specified. Further, the 
friction coefficients have bGcn taken to be constants. 

In the present I'ork they have been treated as functions 
of Re 7 /nolds number, surface roughness and geometr 3 ^. 
Present work brings in a correction to the equation of 
motion formulated b^’- Paidoussis't 48, 61 3 . In view of 
the author, Paidoussis' representation of the norinal drag 
forces on a cylinder subjected to axial flow, though 
accounts for fluid induced damping when fluid is sta- 
tionary’', fails to correctly represent the drag for small 
values of flow velocity, since the drag depends on 
relative velocity between the body and the fluid. 

Since the present work deals with two flows, 
their relative direction for symraetrically supported 
cylinders and individual directions for cylinders with 
asyimnetric end conditions become important. To the 
authors knovrledgc the importance of the direction of the 
flow has not been considered in published literature so 
far . 

Briefly, then, the present work adds several 
more relevant parameters to the problems considered 
before and attempts to generaliEe the parallel flow 
induced vibration problem in addition to emphasizing 
some of the parameters in their true perspective. 



CHAPTER II 


THEORY 


11*1 Introduction 

In this chapter a generalized parallel flow 
induced vibration problem has been formulated. It is 
not being claimed as the most general one since it 
has been restricted to small transverse vibration of 
slender cylindrical elements subjected to steady 
incompressible turbulent flow, para.llel to the coimaon 
axis of the cluster. This formulation synthesizes the 
problems considered by Paidoussis^S 3 Paidoussis 
and Issid into a single problem of vibration of a 
cluster of tubular cj’^lindrical elements carrying the 
fluid. iVhen appropriately specialized it covers almost 
all situations involving parallel flow induced vibra- 
tion investigated so far. 

11,2 Formulation 

For the problem under investigation consider 
a system of identical, straight and uniform slender 
cylindrical tubes (Fig. 2.1), which are assumed to be 
elastic. In practice, the tubes may be held in place 
by end plates or spacer plates at one or both of their 
ends. The entire assembly is surrounded by a rigid 
circular channel through which an incompressible fluid 
floxiTS parallel to the channel axis. It is assumed that 
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the supports do not significantly disturb the flow 
except over a certain length v/hich is small compared 
to the unsupported length of the tubes. Each tube is 
assumed to carry equal flow of another incompressible 
fluid, which may be exchanging heat with the fluid 
surrounding the tubes. This is an idealized model for 
Parallel flow heat exchangers. 

For convenience, the flow in the tubes will 
be termed as inner flov/ while the flow in the channel 
will be referred to as outer flow, Subscripts 1 and 2 
denoting the related quantities respectively. Flow in 
one or both the regions may be two phase flov/, in which 
case the portion of one phase vn.ll be assumed small so 
that the flov/ may be taken as equivalent single phase 
flow with appropriate fluid properties. For the pur- 
pose of expressing fluid forces on cylindrical elements, 
both the flows will be assumed to be perfectly parallel, 
uniform and steady. Further, it is assumed apriori, that 
the vibration amplitude remains small compared to the 
significant length scales of both the flov/s, to the 
extent that the consequential perturbation in the flow 
can be ignored. In the present state of the formulation 
one has. to icealize the flow field since full differen- 
tial analysis of the coupled solid-fluid interaction 
problem is not being attempted, the general form of 
v/hich can be stated formlly as i ^ 
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(a) 

together with 


L '‘lik ( W '^) 9 

(b) (2.1) 

and 


~ D nt /D t 

(c) 


at bounding surfaces 
where 


7^ 1 ^ 1 ^ = kth component of the motion of the 1th 

element 

u^ = ith component of the flow velocity 
p = fluid pressure 

= kinematic viscosity of the fluid 
D/Dt = lulerian Derivative. 

In the above set of equations, (2,1 a) are 
the Navier-Stokes equations for constant property 
fluid, written in tensor notation. Equations (2.1 b) is 
a set of differential equation which must be satisfied 
at all bounding surfaces, the vector operator con- 
sisting of various derivatives to account for flexural 
restoring forces, tangential and normal forces d.ue to 
surrounding fluid and inertia forces arising from the 
motion of the solid. Finally equations (2.1 c) is a 
set of no penetration and no slip conditions, taken 
for impermeable surface and viscous fluid. 
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The set of coupled equations (2.1 ) which 
are generally nonlinear partial differential equa- 
tions, is difficult to solve. For rigid walls and 
parallel flow, the wall pressure fluctuations have 
been studied recently by Schumann '{[6 2 1 by simulating 
equations (2.1 a and c) on a digital computer. but 
the method is computationally uneconomical- Studying 
the motion of the solid coupled with the fluid flow 
by simulation would prove to be much more difficult 
due to the appearance of differentially coupled boun- 
dary conditions (2.1 b, c). 

Returning to our idealized situation, con- 
sider a small transverse motion y of one of the inner 
cylinders about its position of rest. It will be 
assumed that the motion is planar and independent of 
the other spans of the same cylinder so that the 
cylinder need not be treated as a inultispan structure. 
Further, the flow field in two different cells 
(Fig. 2.1) and consequential excitation on cylinders 
surrounding the cells will be assumed uncorrelated so 
that an element can be analjrzed independent of others 
and the entire system as such need not be considered. 
This amounts to assuming that the flow in any of the 
cells is uncorrelated to the motion of the cylinders. 
This assumption greatly simplifies the coupled problem 
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and enables one to consider a single cylinder in e. 
cluster as tiie system excited by the flow around and 
v/itliln. As the spacing between cylinders and signifi- 
cant length scale of the flow increases compared to 
the vibration amplitude, the accuracy.'- of the previous 
assymption can be expected to increase. Angle of 
incidence i and curvature 1/'^ x are assumed to 
be small so that separation in the cross flow is not 
likely to occur and forces on a section of cylinder 
under consideration, exerted by the outer flow may be 
assumed to be the same as those acting on a correspond- 
ing section of a long cylinder having identical cross 
section and inclination to the flow direction. The 
forces and moments acting on a small section cTx. of the 
cylinder are shown in Figure (2.2 a). The assembl;^ is 
assumed to be inclined at an angle 0 to the gra\^ita- 
tional direction. F^^ and F^^ are the resultant longi- 
tudinal and normal fluid forces per unit length, acting 
on the section, while T is the axial tensile force . 

The forces F _ and F„„ which are the result of mean 
px py 

fluid pressure, will be discussed at length later. 

Consistent with the small motion assumption, 
the spatial derivatives of y will be assumed small. 

For the section of the cylinder under consideration, the 
force balance in X direction gives, 
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"XE + Ft. ^ - f 

*bx_ ^ ^ > 'f 


w 


- O (2.2 a) 

where m is the mass of the cylinder per unit length 
and 

^'u =■ '^'1-1+ ^l-x. 

f ^ ■= F^>., 4 Ffvji. 

Forces acting on the section of fluid, contained in 
the aboTe section of cylinder, are shovm in Figure 
(2.2 b). The force balance in x direction can be 
written as, 

.-1. ( K ^ ^ 

where p is the fluid pressure, ^ denotes fluid 
density and A is the cross sectional area of inner 
flow. 

Force balance in y direction, for the 
element of cylinder and fluid, respectively, may be 
written as 


(2.2 b) 


XL 


h -V F u. F b jj 


c 


T 


/ (2.3 a) 



— Irn ( 

^ '6 3 

-3^') + ^ 

and 






A , 

W F A - 

where 

■>v. and 


are inertia 


the shear force on cross section of cylinder and 
f (x, t) is the flow originated excitation force, being 
considered separately from Fjj and assumed independent 
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of y or its derivatives.. Since- the elements have 
been assuiued slender, shear deformation and rotary 
Inertia -will be neglected. Further, if the nmterial 
of the C3rlinder is assumed to follow Kelvin- Voigt type 
relation, the moment balance can be written as 

-1- +,M.i (2.4) 

where is the material damping constant and El the 
flexural rigidity. 




The forces F and F arising from mean 
px py “ 

pressure of outer fluid can be obtained hj integrating 
around the outer surface; however, following the sugges- 
tion of Paidoussis 48 1 j they can be found conveniently 
as follows : consider the elementary section shown in 
Figure (2.2 a) momentarily frozen and immersed in the 
fluid on all sides as shown in Figure (2.3); then, 
there wmll be additional forces A and 
£'^'C p-,a)/3x 3 lx on the cross sectional faces 
of the element, but the resultant of these forces and 
of the forces and V'p'^ dx is knorni and 

equal to the buoyanc^^ force. Accordingly’', 
will be taken at most linear func’tion of x, which 
covers mean pressure distribution due to gravity and 
skin friction drop. The force balance then gives, 
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j ^ A ct C ^ ■>'■*= '=^) 


C Are ^) ^ 

X7 ^ (us'.YN^ Q-waSs’s 

^ ’ D,VfVg#v,C€ 


= "■ ^ 


Cvo\.) 


'll- A I 

X 


Hence , 


\- 


tr 


Itx 


(a) 


(b) 


(2.5) 


It can be observed that the force F is 

py 

equivalent to a tensile forqe of magnitude- (p^ 1) 
acting on the cross section of the cylinder, also 

if = 0; however, for convenience, this 


F = 0 
px 


'd X 

will be invoked later. 


Substituting equations (2.5) Into equa- 


tions (2.2 a) and (2.3 a) 


l==A) ^ A) + li - ° 

(2.6) 

1^ - f M ^ ^ul2. +1 a t t’lA) 

"bx Ot 

(2.7) 

Adding equations (2.6) and (2.2 b) 


('^3x.+ -» hi + F« 'll ' ® 

(2.6b) 

where 


^ A p>^-~ 




32 


Similarly, adding equations (2,7) and (2.3 l) 

1 a k 






(2.7b: 


As will be seen later, if 


N2. 


*&TC, 


OCs’-) 


^ ^ O(^), d *^^07 then 

and rest of the terms in equation 


(2,6 b) are ^ OCs.) ; thus neglectin^g the last term in 
equation (2.6 b) and integrating betv/een the limits 
z and L \diile invoking azial uniformity, 

(t 4 pa)!^-- CT+t=ft)i, 

^ *^ ( 28 ) 

Combining equation (2.8) v/ith equations 
(2,4), ( 2.6 b) and (2.7 b) one gets, 


AJl 4 ?,t) a] (l- X) 

■*■ S ( n»sA'- Ah ^ 

The terms ^ ^ and buo- 

are yet to be determined in ezplicit form. 


(2.9) 


The pressure drop due to gravity and wall 
friction drop can be written as 

( 2 . 10 ) 

where is the bulk velocity of outer flow, 2 is 

the overall friction factor to be chosen according to 
the geometry and Dj^, the hydraulic diR,meter, The 
discrete variable ^ has been introduced to take 
into account the direction of the flow and hence of 
the resulting pressure drop. Thus ■ 
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§ = + 1 if the flov direction is along the 

positive X axis 

= - 1 if the flow direction is along the 
.negative x axis. 

Henceforth, hulk veJ.ocity will always he witten as 
% ^ the variable £ ^ being dropped where unnece- 
ssary. 


The forces and exerted hy the fluid 

shall he expressed in terms of the gross momentum of 
the fluid. The force Fj^. can he thought of as con- 
sisting of two parts; one due to viscous drag, 
and the other due to inertia of the fluid, Thus, 


As proposed hy Lighthill[lO 1, the force Fjj^ is taken 


( 2 . 11 ) 


where M is the virtual tiass of the fluid. 


For the problem under consideration 


F 


TS\ 


^ *"tw1 ' 3^^- 2. / \ , . I > C VI 


''3k\ 

'bt 


( 2 . 12 ) 


where ^ is the added mass coefficient for the outer 
fluid; for inner fluid added mass coefficient has been 
taken equal to unity. This coefficient has been com- 
puted hy Chen and Wamhsganss \_ 67 3 for annular flow and 
hy Chen ^683 for cluster of cylinders. 





For the forces T-^j- and following expre- 


ssions were proposed by Taylor^. 5 J and verified 
Hoerner ^21 3^*03? steady outer flow. 

fvn^ = C I i fj? O 

4 uf -^5 '• 






(2.13) 


where is the wall friction coefficient and 

the form drag coefficient, D is the diameter of the 
cylinder. can be trritten as 

C, ^ ,r r 

^ I = Ub^ -t t) l- 

v/here L is the length of the cylinder. 

The expressions (2.13) can be modified to 
take into account the direction of flow as, 

f'vMi - 4 ?vDUb c tj; s'.wo (a) 


LZ 


= 4 ?^!) '■ 


(b) 


(2.13) 


Equations (2.13) give only mean values in 
terms of the gross momentum of the fluid, since the 
coefficients vary from point to point v/ith the deve- 
loping floW) further, the forces and are 
assumed to be instantaneous without memory. For 
periodic Lairinar Boundary’' Layers the wall shear stress 
is known to lag behind the fluctuating velocity of the 
fluid. Exact nature of the viscous interaction in 
unsteady turbulent flows is not clearly known since 
these flows have not beenstudied in sufficient details 



sun ace 


so far. Yalidit^- of the equations (2.13) for a 
excited by flowing fluid appears questionable f but if 
one is not interested in detailed notion history of 
the surface, it is a simple first approximation and 
offers an attractive alternative to the treatment of 
interaction problem. As will be seen later the force 
plays an important role in the stabilitj’- of 
cylindrical element and hence deserves more detailed 
inve s tiga tion . 


and I 0^ 2 . angle of incidence 5 I - ^ 

4- C ^ 5 will be small and 

approximating sin i by /C^'bJbT^ 

in equations (2.13), while dropping P. from the 
coefficient 

^ V 1 ' fa J fa -h ^ k ^ 

+■ (a) 


I'or the present problem, for small 

1 ^ ' f b 


n 


r Li = 


The second term in the equation (2.14- a) 


(b) 


(2.14-) 


represents linearization for small U-j^ and in this 

form it correctly gives the direction of the normal 

drag force. It can be observed that all other terms 

in the equation (2.l4- a) vanish for U, = 0 and non- 

°2 

zero value of accounts for viscous drag and con- 
sequential damping due to surrounding fluid. This 
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term slightly differs from the linearized term taken 
by Paidoussis fn-S | and in this form accounts for the 
viscous drag due to transverse velocity of cylinder 


relative to the fluid. 

Finally consider the term 


(t t 


For the case uhen the cylinder is free to slide axially 
or is completely free at x = L, recalling that the term 
pA arises from the pressure acting on the sides of 
the cylinder, at x = L '^L ~ ~ ^b -iiere p-^ is 
the base pressure. If the cylinder carrying fluid is 
discharging to the atmosphere at high velocity, the 
reaction force may be taken care of by appropriately 
selecting p-j^. Hence, 

^ ^b') (2.15) 


which is somewhat like a form drag. In case of a 
cylinder free at down stream end the outer flow 
indeed causes a form drag of this type. All these 
cases may therefore be grouped into 

(t-'pA) = if (2.16) 

where is the form drag coefficient to be chosen 
appropriately. On the other hand, if the tota,l length 
of the cylinder remains constant between fixed supp- 
orts, as may be the case when the C3rlinders are held 
between two rigid plates and axial tension is applied,’ 
then, at x = L there mil be an additional compressive 
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force - L > assuming that the weight 

and axial frictional drag are taken by both the 
supports equally. The existence of the second term 
depends on whether the length was fixed while the 
cylinder was horizontal or not. There will also exist 
a compressive force d'ae to radial contraction arising 
from mean differential pressure p. For a thin tube 

this is where . A - bA, and 

’ ' - 2 - 

is . Poisson’s ratio for the material. Finally 
considering externalljq imposed uniform tension T , 


(-r4\=A)| = 




(2.17) 


v., u ' ‘2s'x 

combining equations (2.16) and ( 2 . 17 )? 

(■T+!=.h!L=^'7'-zv)yA -T 




where 


-p { l-F) X) Jh> 

= 0 if the tube end at x = L is free to slide 
axially . 

= 1 if the supports do not allow net axial 
extension. 


(2.18) 


:■ For the problem under consideration equation 
(2.18) can be written ^ ^ ^ ^ 

(T4t=A)l= J[Ci-Z4>)i;A, 

4- 4r ( \-d ') ^ k) (2.1b 

Where depending 

on the particular case being considered. This alter- . 
native is necessary since one of the flows may have 



zero velocity, with ^ - 0 and nonzero force nay eivist 
at end x = L. 

Sutstitution of equations (2.12), (2.14) and 
(2.18 a) into equation (2-9) leads to 




'a'" 






't ( It +'^'=.-v ^ 

a. 


Sv ^ 'N'-j 


'S-t 




c\ + 


1 + [c— ?y ^ t ^ ^ 

4 /Ow] (— yoay '=^ 

^ 


- L 1-^ , . ^ 5^t i /I>K- i ^.] 


o )j 




(2.19 a) 


The following set of general houndary conditions can 
he used for the cylinder, 

E 1 tU ^ V 


■“hTC'^ 




tl 1-1 

~boC'^ 


Co 'll =- ^ 


~?>OC 


tr. 1 r. n 4- 

"b-xl- 


k t I ^ C, = 0 «■ 


k I 


(2. 19 t 


"ax’’-' 


L. 


t 3C, 


where k , k-r and C^, C-. are linear and rotational 
o ’ L o’ L 

spring constants respectively. Depending on the 
values of these all physical boundary conditions can 



Equation (2.19 a) can be v/ritten in the 


form ^ 

T_ : j 

in which the constant term independent of x, j’’ anci t 
has been though of as included in f Cx, t). 

It can be shovni that 


2i-4- - (2.20 3.) 


- ,ui - ^ El ; A 

^ 4- ax 

"^10 =■ LM+ (c,-f,.) ' 3 ^+ a 

‘^xa - f, A uA-i r 0-2^') fh +t] - (,-'-3x7 

7wO-|V-h]- 4_(>-S') ?^-uJ 

All other coefficients *^‘=^'|2 j'^21 22 ’ 30’ 

5 o< 22 j ^ 1 ^ 2’ ^00 ^ considered to be zero. 


II. 3 Dimensional Parameters ; 


In order to bring equation (2.20 a) to a 
diniensionless form following diiner-sionless variables 
are introduced 


^ y '3 I — ( uX l’- t 

L ^ -A(A,^?xVw J "- 

So that equation (2.20 a) reduces to 

il tu 3!!)IL. , 3, (^,^.) 

(To Uo ■> 


(2.21) 



where 

t 


4-0 * ? 


\ =. 
X 


i. t ('f'-'i + P, A+fo *^ 

• 1 


/< 


S I |4*- 

O.' ‘ ^ ^ ^- 


(h,t) to2.- i^CX-i)ji 

P ^ f-x'^/C ^-hf, A +4A) ^ , 

t,,- ^ u 




4 1 ' 


\J I I — 1 

I I I ^1 ^ 

^ fc i 


^ K r - ^ CT-nL-l i - 

to, = r^, ’ 


t,o=iutc,£ + 


-,o- iu..,. 

v^ = D/Ik , e = 'i ^30 , 1^= ^ 

t^o^ uVXu\-^ [(i-2.'rf)n +-r]-Wcos8 44.'j\^^K^'^ 

+ ^ ut£‘i^fc2.}<'-|-?)-£ 0-&) 

w'^ete ji = it/tl , T = T C/*: ■> 

Gb- C, 


the last term in the expression for is 


I 

‘2- 


depending on 


2_v.- 

the particular case under consideration. 

All other coefficients will have zero values 
Equation (2.19'b) can he written in dimensionless 


form 


in X K n = ^ 

•5^3 ^ 


K. 


an 




^0 


-. 2 -. 


^ +w,r\ ^ 11 - k' at (2.22) 

■4-e,3 ' '- 



■where 


, n, = /cti'! 


k'o - Co U /C £ 1} , li. = c /( t I } 


II .4 Response and Stability : 


flo'w, con'ye3ring fluid and si;t jerried -1: : 


The equation of notion of a c;,dui'.fcr in axi.al 

Eli ErlDitrciPY 

force field lias csen obtained ii. article II. '1 . It 
may be t/ritten formally as 

cC j = £ (2-=-3) 

where is a linear differential cperato:.: (equetioii 
2.21). When the force field is generated by the pressure 
fluctuations on the surface of the cyli'idp"'-. the force 
in a given direction is obtained by integrating around 


the periphery of the cj^linder, i.e. 

ah 

=2i= \ K?. h,r) 

-*0 ^ E i ^ 


^ _ lo .4 

Ob g, o'. 




where p is the dif ferenti.al pressure, acting over tht 
surface of the cylinder. In ordrc: to ol ’Gain tnc res- 
ponse due to an arbitrary force field, the c ontinuous 
system is transformed into a discrete one to seeking 


the solution of the type 


Y 


where (t) are the generalized coordinate-s and 
I constitutes complete set of eigenfunc- 

tions of. the free vibration problem of a beam with 


(2.25) 



the same boundary conditions 3 .s the- cylinder without 
any fluid. 


T i C 

J r. *- ' ' 




-iC ua ti on (2.21) the n b o c ones 

T. EE . 

T 1-0 3^0 ^ 

In exnanded forn this becomes / . 

E Pr i t„ q., 

+ tio (2.26) 

where dots denote dif ferantiation i/ith respect to 
and primes those with respect to 'K . 

')k,^ are the bcs.m eigenvalues corresponding to 
eigenfunctions satisfying, 

For the boundary conditions considered 
he rein \P- . ^ is a self-adjoint set such that^ 


(2.27 a) 


0 

and 


( <p^<^ cl^ = Sys 


y 

0 


(2.27 b) 


(2.27 c) 


where A is the Kronecker delta, 

YS 

tfultiplying equations (2.26) b 3 ^ a i 

O 

s = 1, 2, ... and integrating between the limits 0 
and 1 vields 


z 

T 




(2.28) 



I 


where 


St 


- t 

W -I K.J 


az. 


Qs - ^ f, c^ T) 


0 


Trxxncn.ting the sucimation at r = H; equatioi 
(2.28) may be written in marrix form 


H 5, + K 1 = Cl 


(2.29) 


n can be chosen to be large enough to determine 
the response -v.dthin a desired accuracy, latrix M is 
diagonal but C and K are asymmetric matrices due to 
the presence of Coriolis force terms. Equations ( 2 . 29 ) 
represent a set of coupled ordinary linear differential 
equations idiich is difficult to decouple unless matrix C 
is a linear combination of matrices M and K, which can 
occur only in particular cases. For a general case, equa- 
tions (2.29) can be reduced tc a system of first order 
differential equations as follows ^ 

■i' 

Pre -multiplying equation (2.29) by IT ' 


^ C ^ + K K 


(2.3c; 


let 


A. 

Z d= 


5- 




From equation ( 2 . 30 ) one gets 


i=- & 


C 




(2.3^ ) 



where 


K ^ 


0 


I 


-m-'k ! -m‘'c 

{ 




J £'2-1' \ 


o 


H*Q, 




where £ 0 ] is the null matrix and £ I 3 is the unity 
matrix. Equation (2.31) represents a 2n dimensional 
linear time -invariant dynamical system with system matrix 
B. Matrix B is generally asymmetric and nonsingular. 
Premultiplying both sides of vector equation (2.31) by 


-1 

B , which is presumed to existy 

I) I - z -h df 

where 

D = B 

For the rector equation (2.32), if one seeks solution 
of the type 

Z = e u 


(2.32) 


the eigenvalue problem 


[Dm -I 3 u = [o] 

is obtained. The adjoint problem corresponding to 
equation (2.33) being, 


(2.33) 


[ Dh-l]i = [Oj 


(2.34) 


In general, the eigenvalues of equations 
(2.30) are distinct and it can be assumed that 2n 
eigenvalues of equations (2.33) and (2.34) and asso- 
ciated modal matrices A and N can be obtained such that 



Kl'^D A = j 


and I are diagnoal r^tr^ 

Reverting now to the non-homogeneous ecaaolons 
(2.32) and letting 

5 • * 5 <=.3» 

equation (2,32) hecomes 

DAt, = AC, r DF 

T 

premultiplying "by N , 


jt = + !> 

where 

f 



Since matrices <1 and L are diagonal the system 
is decoupled. If generalized excitation force is 

known in deterministic form, the response to the force 


field may be obtained as 

MlZ 7 /‘l(r-h)x 




o 


cp- ^ ^ 


(2.37) 


and using equations (2.25), (2.31) and (2,35) can be 
transformed to dimensionless response T\C^ 

The first term on right hand side of equation (2.3?) 
represents response due to initial conditions or free 
response alone . If the system has snail perturbations 
from its quiescent state, eigenvalues ^LLi will 
determine whether the response remains bounded. 
Eigenvalues JU-i will in general be complex and their 
real part determines the stability of the system. 
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iifiien the force field is random tut stationary 
and ergodic, aiierein lies the chief interest of this 
nork, the response statistics are ottai.ined as follows. 


The space tine correlation function of the 
response is defined as 


Jo) = r ), + s,,r+T,) > 


— I I'm 

T-m «= 


which, in view of equations ( 2 . 25 ) (2.30) maj’’ te 

written as ^ 

-r 




— llVYl _L. 

-X-j^osT 

o sj 

where n is the order of vector equation (2-30). 


Eqiiation (2.38) can te transformed using Correia. tion 
The ore e: to yield 

■ * 4,, •- J (2.39) 




(r^ 


where ZC03> is the Fourier Transform of z.. v e ^ 

I'Jext, taking Fourier Transform of equations (2.35) 

and (2.36), respectively, 
zv. 






and 


ibd J Y — L^x ^ Z 


V ^x ^ ^ 


Z, C^) - 



(2.40) 



where the Fourier Transforsi of 6-s Cf ^ 

•r\^- 

Also it is convenient to define 

H-d' C^O ~ i u) Ty ~ Ly ^ 

Substitution of ecuations and (S.'m-I) 


-} 


V£:, I 


in equation (2.39) gives. 


2-y<.. 


* ' t, • * ^ ^ Q L’ ^ 


“J" — :$> 0*3 T| 


— oo 


Defining generalized excitation cross spectral density, 


s.h-) - 1; 




T 


‘f 5^*^) (2.^3) 


equation (2.42) may be rewritten as, 

V”. (:■?, ?o.'T.j = -4. f_ t/su,)L '■hv 4->' 

I J : • ^ '»' r ^ 

( C-*^) ^ ys ‘^a) c(‘^ (2.44) 

—=io 

From which mean square response nia.y be obtained b 3 r 
setting ■ 

Thus, 


■n 






zrt 


> 


t j 


2.’vi 

Atv a_,-3 

' (2.45) 


Generalised excitation cross spectral den 
sitjr needs to be examined in detail 


! i ry\ J— C~'^') ^ s 

r-^ on T f \ 


•vH 


i 


I T P 

wo ^ X 


!;>>< _1- ^ ^ V^y_ntk 



where P 




T 

= n E‘ so that 

- i ' 


- liW 

k,JL 


Y T-*'«3C T 


U ft K, \ _ I. 

AXM 


i- -2_ 


In equation (2.46) the quantity in square hrackets is 
the cross spectral density of dimensionless excitation 
on the cylinder. In order to study it further, consider 

the space time correlation function \iXf , given by 

T 

o 


, 1 Tit ITT 

, \\y^ _L 1) L ( (( 


, ~ \ e t)K'^eu,T+ro) ccsy,cc£B.j^M^i 

O o<^ 

Interchanging the order of integration, which is per- 
inissible since integrals are assumed to converge 


r 


■“o^T S P' ' 


^2- C t^r XTt _ ,. 

\lJj Ck, - ^j=- ( f 

^ o o o 

t.tr 

"— C( "'vlj c ^ 'r 6 ^ T o')dOs0,C0sB2<^®’^'^^2- (2.47/ 

£> 0 


Using Correlation Theorem, equation (2.47) can be 
transformed as 


I L eo Co , 

\ 'iw — ^ —Ini') v (f T' tb) duj 

-p^oo itTT ) ^ 

CJ-O. 


C ^IsXj. 


iyvi 


-U£ r r 

C 


i<2>ro 






Clianging the order of integration 


\ \ YY\ \ 

"f "4 552) 2 Tv ’ 


o<j 



I 




— ow 


riT 


= ^ [-tfii ( ( '■■ 

'*-£'■1'^ J 't. 


O 0 


rr-n. .6,-^) 'pC-^^, 

xr\ i I •» I y / . 2. / 


Sx wy 


OOsOj Co? 9’2 _cI ^ 2 . \ J.jj 


Comparing both side 

\ \VY\ _1 

T-^ oo T 




ITT ZtT 


' T— a ca 
~ 0 ^ 






%Tt 2.Tr 


pT /ij) CO^ &^foi 


n r 




^ -5 


0 c 


^w) CCS 6 ( C6? 

(2.48) 


■where P pp C"^ v tx / cross spectral 

density of the differential pressure y In the experi- 
mental results available so far, this quantity has been 
studied extensively. 


(2.46) 

■P ( 

ITj ,^tT' 

U 

c o 


Substituting equation (2.48) in equation 


! I 
( 


pp 

4 1 (> 


L r ^ 




P, 

, »^ + k 

•4 g *- T X. L — 

k,-i 


^ T & 

■> ' S X , ^ i ,, 


(jJ ) 

« > 


I \.) 

Cos 


n 


t=. DO A 




■5- 




(2.49) 


k , 1 



where 


I { 


^ ^ ^ 1= C?, , <^. 


ki 

and 


O <5 


( 2 . 50 a) 


2F iTT 


' O O 


Thus, for the case i/hen e::c5-tation is due to "•.'a-ll 
pressure fluctuations, the space tins correlation function 
or spectral density of the response can he obtained if the 
vra,ll pressure correlation function or cross spectral 
density is known, along with various coefficients used 
in expressing fluid forces on the cylinder. The latter 


must bo specified in order to obtain either eigenvalues 
or response. The response obtained here is a linear 
mapping of random wall pressure field into idealized 
twfo dimensional response field. In the excitation 
considered here, the fluctuations in the wall shear 
stresses have been discounted since their component 
in the direction of motion of the cylinder, in an^^ 
case, be negligible. In actual practice the motion 
will not be planar due to Coriolis forces and some 
torsional deformation will occur. In this f ormulaticn, 
other tji^pes of excitations, e.g. structure -borne vibra- 
tion, can be accomnodated easily. In this work, however, 
chief interest lies in stud 3 ring suability and responso 
due to turbulence originated excitation alone . Other 
types of excitation like structure-borne vibration, the 
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parametric one arising from pulsations in the flow 
rate and consequential vail pressure, are characteris- 
tics of each particular situation and therefore offer 
little motivation for the study. 
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CH.4PTER III 


THE FLUID IIJTERiiCTIOII PAEAIETERS 

111.1 Introduction ; 

In order to solve the equation of motion (2.20) 
developed in Chapter II , for obtaining free vibration 
frequencies and associated damping or response to arbitrary 
forcing function, one needs the fluid interaction coeffi- 
cients G^, C^, and excitation f^ (or its statis- 

tical description if the excitation is random). In this 
chapter the fluid interaction coefficients C^, and 9C 
have been studied and expressions have been derived for 
them from the available experimental and theoretical 
results. 

111. 2 Normal Drag Coefficient ; 

The Normal Drag Coefficient C^ has been introduced 
to take into account the viscous drag vdien the outer flow 
is negligibly small or zero. For a general ease considered 
here, this factor will depend on the geometry of the 
flow passage and Reynolds number. For a general geometr3?' 
this factor has not been studied in detail either theore- 
tically or experimentally. In the absence of direct 
measurements one can obtain the average value of this 
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coefficient from the measurements of viscous damping due 
to stationarjr fluid or at small flow rates. Such 
measurements have heen made "by Chen and l»ambsganss [*67 ‘j 
Pavlica and I'arshall ’ 23} .32| and McCormic and Kippley 
Yp9 J • Here an attempt has been made to derive the 
dependence of this coefficient on fluid flow induced 
damping measurements at small or zero flow rates. From 
the analysis, the coefficient C^-, has been evaluated for 
annular flow geometry. For simplicity, the coefficient 
will be assumed to be independent of Fceynolds number. 
With these assumptions consider equation (2,21) with 

fo = ^1 =TT=r = 'y =0andS^^ = 1, 




.5 




where the coefficients t’ are evaluated with u^ 




Su= I- 


-V-CJ, 

Let us assume a solution of the type 

'^(Ih) = (3.2) 

Substitution of equation (3.2) in equation (3.1) leads to 
« p"® + F ® + i.wt';’, F i- F 1„ F i F = O (3.3) 

For the case lAon velocity U2 — 0, let 

f =. R -t • • ■ (a) 


60^ -t + 603,4- • ■ • 

A " 

. # 


(b) 

(c) 


(3.^) 


where is the damping factor associated iz-ith tne first 


mode. Substituting equations (3.^ a,b) in equstion (3.3) 

o 

equating the coefficients of the terms containing u^ 








locco, F, +5 


or 


where 


B 




o 

Z 


D. Bo =-=^ 


(5.5) 


I? 


1^0 - / O'** 

Equating the terms containing u^ , 

* r" ^ ^ t W O Vy 

^ P| ■+i,ics p 

F3+2U(,^i't'oit " i 
F„-i F, 


i.e 


=-[-i«u,F„l-* + e2<0o“it,r-h“.':+f'"<£ (3.6) 

- i ito.c fi p - (2 tX ^^ 0+ ps) '"/J 

Since must satisfy the sane boundaiy condi- 
tions as those for 4 > in equation (4.2) 

i, 2, . . - 

From this we obtain the first mode complex frequencj", 
using equation (3.5) 

The damping factor associated with the first mode at 
Ug = 0 and for the case c<= 0 is now obtained as 




ce f 






for all boundary conditions (3.7) 


Multiplying equation ( 3 . 6 ) by and integracin 
with respect to ^ between the liraits 0 and 1 leads to 

iC 

if F (o) = F (1 ) = O 5 then for the case oC = 0 , the 
contribution to the first mode complex frequency is 
given by 

CO ^ 

g 3^- i t.c e f/^ 


which gives the contribution to the damping factor asso- 


ciated with the first mode 






If F(0) = 0 but F(1) 0, then 


(3.8) 



From equation (3.9) it can be observed that the Coriolis 
force term can promote damping only when displacement at 
one end is nonzero. 


In order to find the dimensionless parameter C 
and consequently viscous drag coefficient the first 
mode damping factor measurements by Chen and Wambsganss 
[^672 were used. These measurements were made by 
a) bandvfidth method, b) magnification method ’/itn cons- 
tant force and c) magnification method I'.dth constant 
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:.isplacement. The dependence of the first mode damping 
factor on the velocity of flow was measured hy applying 
external electromagnetic excitation to a cylindrical beam 
with clamped-clamped and pinned-pinned ends. The method 
of applying excitation did not disturb the flow and the 
displacement of the beam was observed by optical means 
for three different confinement ratios ( n = 0.33, 0,5 
and 1.0). For the annular flow geometry these experimen- 
tal results were relied upon. The curves of the first 
mode damping factor against the velocity of the flow x^ere 
approximated by a quadratic polynomial for the inteilha.1 

^2 ~ ^ 6.55 and an expression for was obtained as 

C X) 2 

= £ (^t20-SC73cl,2r^’3l*t36o!d,2_4- 

where d^^ = radius ratio = h/(1+h) 

^2 ~ dimensionless kinematic viscosity of the 

fluid 

= kinematic viscosity > /3£l) 

V EX 

The above expression is applicable to a single 
rod supported in a circular chs-nnel, for the range of 
confinement ratio 1.0^ 0.33. 

For an isolated cylinder (d^2 
computed from cross flow normal drag ■coefficient, which 
for low Fieynolds numbers, can be approximated as 

^ {O ' o/Re 



where Re 
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Reynolds number based on the diameter of 
the cylinder 

Cjp = Drag coefficient = Drag Force /I pu^£;;L 

from which 


r rv -(I 


so that 


c 


0 .- 


V 


z-j-2.ev. 


where i>. 


dimensionless kinematic viscosity of the fluid 


U 


^2. A 
E X 


Since Cd has been assumed to be independent of Reynolds 
number, following form of functional relationship has 
been assumed 


u 


or constant 


In the above relationship the constant on the 
right hand side was assumed to be 10.0 for an isolated 
cylinder. 


Since only a space averaged value of the coeffi- 
cient Cp is required to account for the normal form drag, 
these expressions were thought to be adequate for further 


use . 


™*3 Qkin Friction Goefflr.ipnt. 

Overall friction factor /\ has been measured 
by numerous researchers for various channel geometries. 
Direct measurement of the skin friction coefficient is 
difficult. Both coefficients vary \jitli Rejmolds number but 

r-' 

for Reynolds numbers greater than 10^ the variation is small 
The measurements of friction factor for annular flow by 
Lawn and Elliott |_70J and Rehme indicate that the 

overall friction factor is greater than that for a 
circular pipe of equal hydraulic diameter by upto 5?^ 
(depending on the radius ratio d^2^’ Reynolds 

numbers. The skin friction coefficients for inner and 
outer wall, and respectivelv, can be found if 
zero shear stress radius is known. Since the flow is 
radially asymmetric the radius of maximum shear stress 
and the one of maximum velocity differ, the difference 
being more for smaller radius ratios, F^ehme *^71J collected 
the experimental data of several other researchers for 
annular flow and found that these radii correlate as 




where = maximum velocity radius 

d = zero shear stress radius 
o 

d^ = diameter of outer cylinder 
d^ = diameter of inner cylinder 
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The skin friction coefficients can he estimated hy usin^ 
the laminar flow ratio 


using equivalent friction coefficient for the annulas 

r - 

^ a, -t a^ 


( 3 . 11 ) 


one gets 

Cf^ _ I — / 4 ^ 

I — 

which gives the inner wall friction coefficient C^, 


For computational purposes, the overall friction 
factor ^ for confined flows was first found from 
Colehrook-White Relationship 

j=- IV'f-0-87k(^ + (3. 

where k = equivalent surface roughness = _J — t — 




( 3 . 12 ) 


d cA 2_ 


for an annulas 


for a cluster of rods. 


d^ denoting rod diameter and d2 the channel diameter. 

The friction factor ')'» obtairied by equation 
(3.12) was increased by for annular flow geometry 


irrespective of radius ratio, since exact dependence of 
^ on radius ratio is not yet clearly established. 


m Added Mass Coe ff 1016111 : 


The added mass for a body moving in a fluid 
arises from the hydrodynamic force required for accele- 
rating the fluid replaced by the body. It is somewhat 
like a form drag in inviscid fluid. Generally, the 
added mass of a cylindrical rod is assumed to be eq.ual to 
the mass of fluid replaced by the rod. This is true for 
a rod submerged in an infinite fluid. However, in a 
cluster of rods, the added mass is affected by adjacent 
rods and channel vra.ll. The added mass increases v;ith the 
increase in confinement. Chen and Wambsganss \l673 have 
computed the added mass coefficient ^ single rod 

in a channel. Chen C.683 has shown that the added mass 
for a cluster of circular cylinders is a second order 
tensor. The added mass coefficients can be obtained in 
the form of a series solution. In general, the self added 
mass coefficient for the central rod in a cluster is much 
larger compared to the one for the same rod in an annulas 
of equal confinement ratio h. The results of Chen and 
Wambsganss show that the added mass coefficient is a 
strong function of confinement ratio h but weak func- 
tion of sound velocity and frequency. 

In this v/ork, for the computational purposes, 
the added mass coefficient vras computed using the expre- 
ssions given by Chen and Wambsganss |^67 J for a singls 
rod in a circular channel. For an isolated GVj.maer 
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this coefficient was taken to he unity xdiile for a cluster 
of cylinders this coefficient was either assumed or talien 
to he the same as the one for an annulas of equal confine- 
ment ratio, since the explicit expressions for the added 
mass coefficients in terms of series solutions have not 
heen given hy Chen 

Since the response of- the cylinder due to the 
flow originated excitation was not computed, the statis- 
tical model for the random pressure field at the wall was 
not necessarj^. Next chapter briefly outlines the computa- 
tional m^ethod used to study the stability of the cylinders 
using the fluid interaction parameters studied in this 
chapter. 




CMPIER vr 


COl-iPUTATIomL APmOACH 


This chapter briefly deals with the computational 
method used to obtain the numerical results. 

IV. 1 Characteristic Equations and 
Eigenfunctions ; 


Solution of the governing equation (2.23) using 
Galerkin hfethod requires the set of eigenfunction to be 
used as coordinate functions. These were computed from 
the ’’abbreviated” eigenvalue problem, 


w 


- 


h- 


(4.1) 


Subject to the boundarj^ conditions (eou-stion 2.22) 


-I- K.4 = <p =0 at f = 

cp'" + k, cf> = (p "c 


D 


(4.2) 


Eigenvalue problem defined by equation (4.1) o.nd (4.2) 
is self adjoint and yields a set of eigenfunctions t 
satisfying 


and 


(45 ^ ^ 




i 


o 


r “s 


1 qp ' -of ^ 

(4.3) 

■ V, , 


M y^YS 

(4.4) 


For the actual Gomputation solution of equation (^.■’) 
is taken as 



67 


■+ CAG'sh 


-4V^'£h AJ (4.5) 
Substituting equation (4.5) into equations (4.2) one 
obtains a system of linear equations 




4' I > 

|K,<is>+X 4, 5»fi 4 c5sO 

Si..X 1 "^X"CosKX ‘-i'X^SmK'X 


k 


k, 


kjCo^pV |"-!XQ>^X ) I 

. ! *4" G ~ St^WX 


^1 


c 


L’^J 


Q 


0 

o 


6^. 


(4.6) 


For nontrivial solution, setting the determinant of the 
coefficient matrix in equation (4.6) to zero yields the 
characteristic equation "which may be "written as 


A 


,5 k/, K,, k/A) =o fr-7) 

The roots of the transcendental equation (4.7) 


can be found for any desired number of eigen"values . Once 
the eigenvalues are found (n = 1, 2, 3) can be found 
in terms of C^. In order to have definite n"umerical values 
for C^, the orthonormality condition for the eigenfunctions 


can be imposed. Thus, 


=l 


^ = b-' 
^ ti 


(4.8) 


which simplifies further 'computation as well. For 
various standard boundary conditions equation (4.7) 
simplifies as follows : 
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-Pinne d End Conditions : { K - 1 < ■=. • K K o) 

X I J Q I J 

Characteristic equation A(>v) -^iA>. (a) 

eigenfunctions ^ ^ . .(b) 

2 . Glamned-Clamned End Conditions ; - i^d 

Characteristic equation A (C^) c |- (^^ > C^l., > = O (a) 


and eigenfunctions (unnornialized) , 

■+'^a. O js'X,^ -+C3Si.vvUX^ -4-Q/^C^k'X^ (h) 

Co5 '"Ojslf' ^ 


(4.10) 


whe re 

^ i ^2." ~ i i^b' ~^\ i 

3 • Clamned-Free End Conditions ; 


^ u/i 3v — tV K3v 
k ' = OQ K, 


k| 


Characteristic equation A(>.) r I + C?iS > Q>sk?^ ^ 0 (a) 

and eigenfunctions (4,11 

4(^)n C, 5,-^^ ^ (t) 

Oq* ~ I ^ (-2.^ > ^3~ / 1 ~ CcsTv. -+ G3S.la 3\ 

4, Partially Fixed End Conditions ; ^V4 — 

When ends are fixed in such a mariner that the 
transverse displacement is completely restricted hut slight 
rotation is permissihle, the situation commonly found in 
practical assemblies like shell and tube heat exchangers, 
fuel rod bundles etc., this is appropriate boundary 
condition.' For this case, characteristic equation 





o 


(a) 


and eigenfunctions 
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~ CL/ Swv ^ H" C 3 , si v\W^<^ -4 - Cl|. C‘S ^ (t) (4 

where C^- I , ^ = 

“1 C3 = 2V'<; -Ci 

when ^ CO , equations (4.12) simplify to 

corresponding equations (4,10), as they should. 

5 • One End Partially Fixed, other Free i ( ^ \ - ^ 

This boundary condition can be used for a 
cylinder free at one end and fixed at the other end in such 
a way that the transverse displacement is prevented while 
allowing some rotation. 

For this case characteristic equation is given by 

^ 7^ (ViViT^CdsKT^ = o (a) 
and eigenfunctions by ( 1 ^ 

4'(^) = C, Sw'X f + Cids ^ e + Cj £ ^ (t, ) 

where r , ^ , 

^^=1, C2=--i 

-4” (T^sV* “t 

d - 4 

Equations (4.13) simplify to corresponding equations 
(4.11) when K^- co 5 as they should. 


and 


Cj - C^Qis'X 
C^= 2.>'/i 
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In this work only five end conditions of the 
above type were studied. The eigenvalues were co'V'utecl 
using a combination of internal halving and linear inter- 
polation method since eigenvalues are well sepa.rated. 

The coefficients ^ were normalized according tc 
equation (k .8 ) . 

IV , 2 Galerkin Matrices ; 

For computing Inertia, stiffness and damping 
matrices appearing in equation (2.28) one needs to compute 
the integrals ^ ^ ^ . 

the last one due to appearance of a linear term in 
in the expression for the coefficient t^Q in equation 
(2.21). These were obtained using exact closed form 
formulae vvO-thout resorting to numerical quadrature. In 
order to retain the desired precision, the eigenvalues 
i = 1, 2, ... and coefficients C^, r = 1, 2, ... 
were computed using double precision arithmetic. In 
order to expedite and thus economise the computation, 
following relations for boundary conditions of type 1, 2 
and h were used. 
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-0 a >' ^ 2 ^ - j ' 0 ' ^ J ' 

= G^- + 2 Glj 

TV System Matrix and Eigenvalues ; 


c) (i^.1 


The fluid parameters such as coefficients Cj 
C^5 X % were computed as described in chapter III 
and dynamic matrix B appearing in equation (2.31) ^ms 
computed using these. It was then inverted using 
Gauss-Jordan technique so that the eigenvalues to be 
computed would appear in the desired sequence. For com- 
puting the complex eigenvalues of the I'fetrix D two 
different methods were used. VJhen both real and imaginarj^ 
parts of all the eigenvalues were nonaero a method which 
is a combination of Mises pov/er method and Hotelling’s 
deflation technique 72, 73l, appropriately specialized 
f or complex eigenvalues and eigenvectors, was used. Fne 
matrix deflation technique reqiTires construction of sweep- 
ing matrix. It was computed using Gauss Jordan technique 
for matrix inversion, to retain the significance in 
desired eigenvalues and eigenvoctor. .agnitude oi the 
eigenvector was controlled during ever3' iteration so tnst 
the quantities did not grow rapidly and precision vas 
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retaiiiGd. Since only first three modes mqt& corisiderec- 
power method proved fast? further ’•’hile computing 
eigenvalues .one obtains the modal mstrices for the S 3 ''stem 
and adjoint system as well. This enables the computation 
of response due to arbitrary force field by modal analysis# 

i#ien an eigenvalue was pui‘e real or irraginary 
the eigenvalues were computed by a different method. The 
computational procedure switched to this method as soon 
as one of the eigenvalues approached close to real or 
Imaginary axis. In this method eigenvalues were 
computed by QR double step method and eigenvectors by 
inverse iteration. The matrix was first scaled by a 
sequence of similarity transformations so that the abso- 
lute sums of corresponding rows and columns are roughl:/ 
equal. The scaled matrix was then normalized so that 
Euclidean norm equals unity. The main part of the compu- 
tational procedure consisted of the reduction of the 
matrix to upper-rlessenberg form by means of sirailarity 
transformations using Householder’s method. Then the 
QR double step iteration was performed on the ilessenberg 
matrix until all elements of the subdiagonal converge 
to zero or in modulus less than 2”^'^ 1| Hilg* 7he 
eigenvalues were then extracted from this reduced form. 

For computing the eigenvector inverse iteration v/as 
performed on the upper Hossenberg matrix until the 



absolute value of the largest component of the right 
hand side vector was greater than 2^'^/(100 n) where n 
was the order of the matrix. After this bound was 
achieved one more step was performed but at each step if 
the computed residuals were groatcr in absolute Vctluc 
than those of the prei/ious step then the vector of the 
previous step was accepted as the computed eigenvector. 

The principal aim of the computational program 
was to compute the response due to turbulent flow induced 
excitation, therefore Galerkin method was used in spite 
of other simpler and faster methods to compute the 
eigenvalues alone . The stabllit3^ study was thouglit of 
as a part of the general programme for the study of flow 
induced vibration of cylindrical structural elements. 



CHAPTER V 


RESULTS Al-D DISCUSSIOU 
V.1 Introduction ; 

This chapter deals irith the numerical results 
obtained for some of the cases of practical interest. 
Dimensionless complex frequencies of the lowest three 
modes have been found as functions of the velocity of 
internal or external flow. The mechanism underlying the 
onset of various types of instabilities has been discussed 
from the analytical point of view. The significance of 
Coriolis forces has been examined. 

In all the figures of tiiis chapter, square roots 
of real and imaginary frequencies have been plotted on a 
linear scale and only magnitudes of imaginary part of the 
complex frequencies have been considered, so that frequen- 
cies of all three modes can be shown in a compa.ct figure 
for convenience of scaling the results. In cases where 
parameters have been mentioned, they were held constant 
irrespective of varying Reynolds number, lor all these 
studies dimensionless viscosity eas been tahen t.o re 1' 
and densities of inner and outer fluid to be equal, 
qjher ever necessary, the loci of nhe roots \hiich lie on 
axis or overlap each other are shown scheciatically among 
the axis, for the sake of clarity. 
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Figure (5.1) shows the loci of the dimensionless 
complex frequencies of a pinned-pinned cylinder with exter- 
nal velocity, U2) as parameter. The other parameters are 
?> = o.i,e.c^=S= x-t , k = / = r = Tl-C =0 
(see reference 4-8, figure 3). it can he seen that small 
flow velocities act to damp free oscillations of the cylinder 
As the flow velocity increases the system may hecome unstable 
in all the three modes. The locus of the first mode bifur- 
cates on the He (w ) axis, one branch receeding from the 
origin and other approaching and eventually crossing it at 
u^ 3.14-; this evidently indicates buckling instability - 
(ch = 0). The instability associated with the second mode, 
on the other hand, which occurs at 5 * 99 j is oscilla.tory 

As the velocity is increased further, the locus of the 
second mode bifurcates on the Re (^ ) axis at U2^=^-6.27. At 
slightly higher flow velocity, U2<=^ 6.33? the positive 
branches of first and second mode loci coalesce and leave 
the Re ( ui ) axis at a point where Re (^ )> 0, indicating 

the onset of coupled mode flutter. The other branches of 
these modes behave similarly, leaving the axis where 
Re ( <0. Since the frequencies were computed for 

U2 4 9 the third-mode buckling observed in reference [ 4-8], 
(which occured at ^ 2 ^ 9*4-6) does not appear in Fig. (5*1 )• 
The second mode flutter, which sets in at U2^5.99 is 
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notably absent in Fig.(3)of reference C ^8 3 , according to 
which the system loses stability in the second node first by 
buckling and then by coupled— node flutter as opposed to 
ordinary flutter . These differences are due to correct 
manner of accounting for tlie frictional forces in the pre- 
sent investigation. 

When fluid parameters are allowed to vary, 
similar results were obtained, as shown in Figure(5»2)j 
except that the critical velocities are different, evidently 
duo to the inclusion of parameter c, which has been taken 
equal to zero for Figure (5*1). 

In order to study the effect of internal flow the 
cylinder corresponding to FigureC 5.1 )has been considered as 
a tube with internal flow with velocity u^ = 0 . 5 . The 
complex frequencies of the lowest three modes are shown in 
Figure (5.3). The fluid parameters were held constant. It 
can be observed that the dynarracal behaviour is similar, but 
the critical velocities are now smaller. The results, when 
the fluid param.eters are allovred to varj?' with the flow 
velocity,^ are show in Figure (5*^). From Figures (5*3) and 
(5.^)? it is obvious that .the internal flow destabilizes 
the cylinder. Parameters of Figures (5.5) and ( 5 . 6 ) 
correspond to those of Figures (5*3) and (5.^) respectively, 
except that the direction of the internal flow is 
reversed. The character of the frequency diagrams is 



similar but the critical flow velocities for oscillaroif 
type of instabilities are higher. The internal flcv in 
direction opposite to the outer flow, thus, appears to 
stabilize the cylinder; critical velocities for buckling 
however remain unchanged. 

The dynamic behaviour discussed above is not 
unique. The complex frequencies for a pinned-pinned cylin- 
ders with u^ as parameter and p = o-4-8 , 

are shown in Pigure(5.7)* In this case too, it can be seen 
that small flow velocities tend’ to damp the free motions 
of the system. As the velocity is gradually increased, 
however, the loci of first mode frequency bifurcate on the 
Re (^ ) axis, one of the branches crossing origin at 

(as for the case when ^ = 0,1), indicating the 
onset of buckling instability. As the velocity is further 
increased the second mode gets less damped and the loci 
eventually cross the Re (o) ) axis, indicating loss of 
stability lij ordinary flutter, at 6 . 27 - Almost 

simultaneouslA^, one of the branches of the first mode 
crosses the origin, two branches coalesce and leave the 
Re (63) axis at a point where Re (^)< 0. Thus the flow 
- velocity wiiich corresponds to loss of stability in the 
second mode is approximately the one for which the first 
mode regains stability. The coupled mode flutter occur- 
ing for the case when ^ = ,0.1 is now absent. In 
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reference C ^8 3 this case has been studied v/ith slightly 
lower value of ^ such that ^ = 0.25* In Figure 4 

of reference £ 481 , ordinary flutter of second mode does 
not appear and system loses stability in common node 
flutter, the role of second mode being replaced by second 
branch of the first mode loci. This difference can be 
attributed both to higher value of in this study and 

to correction in the representation of the frictional 
forces . 

Results obtained for the same cylinder when fluid 
parameters were varied, are sha-zn in Figure(5*8). Internal 
fluid has been assumed to be stationary. It can be seen 
that the general character of the frequency diagram does 
not change except that the critical velocity for first mode 
buckling is slightly greater than 3*14. The velocity at 
which the second mode flutter sets in does not significantly 
change. Parameters in Figures (5*9) and (5.10) correspond 
to those of Figures (5.?) and (5.8) with the difference that 
the internal velocity is 0.5 in the same direction as the 
external velocity. The critical velocity for the buckling 
instability is now reduced while the critical velocity for 
oscillatory instability remains approximately the same as 
in Figures (5.7) and (5.8). Thus internal flow destabili- 
zes the cylinder. Figures (5*11) and (5.12) correspond to 
Figures ( 5 . 9 ) and (5.10) but the direction of the Internal 
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flow is now reversed. It is seen that the character of 
the dynamical behaviour is similar but critical velocities 
are lower compared to the case when internal flow is absent 
and the reversal in the direction of the internal flow 
further destabilizes the cylinder. 


V.3 Pinned-Pinned Cylinder in Cluster ; 

Next, the dynamical behaviour of a pinned-pinned 
cylinder with increasing flows in a cluster can also be 
examined. If the internal flow is absent, this case can be 
considered pertinent to the fuel bundle assembly in a 
nuclear reactor. If internal flow is present this case 
represents the kind of situation in a parallel flow heat 
exchanger. 

A complex frequency diagram is shown in Figure 
( 5 . 13), ■■•corresponding to the parameters (3 = 0.1, G = 0.25, 
9C-4-, ^T- I ^ >*= ■p's- TT-sO , It can be noted that 

the dynamical behaviour is quite similar to tiic.t of an 
isolated cylinder (refer to figures 5*1 “ 5*c). Lowexor, 
the flow induced damping is smaller in this case ana the 
instabilities occur at much lower values of flow velocities. 
Accordingly, the onset of coupled mode flucter appears 
tmce with buckling instability interposed in between. 

The low value of critical velocity for buckling in first 
mode, u^ = 1 .57 is especially noteworthy. The onset of 
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second-mode flutter at ix^-^ 3.095- is absent in Figure (5) 
of reference C^8 3 , presumably because of different manner 
of representation of frictional forces. Similar argument 
can be given for the second appearance of the common mode 
flutter. The critical velocities (1.57, 3.1, 3.1^+, 4.715) 
for various instabilities shov/ that tlie close spacing 
severel];- destabilises the cylinder in a cluster. This 
comes from the effect of increased added mass as will be 
discussed later. Figures (5.1 4) and (5*15) show the complex 
frequencies when the external and internal flow is in the 
same and opposite directions respectively. Internal flow- 
velocity u^ is held constant and external flo-w velocitj^ U2 
has been varied. Internal flo-w is seen to reduce the cri- 

I 

ti(gal velocity for buckling in first mode (1.55 as against 
1.57). Irrespective of the direction, however, the critical 
velocities for second mode flutter and the first and second 
common-mode flutter are no\j lovrer. Th-as internal flow 
further destabilizes the cj^-linder su.bjected to external 
flo-w. This res-ult is of practical importance since, in 
order to improve the heat transfer in parallel flow heat 
exchanger, higher flow velocities are employed. The effect 
of the added mass on buckling instabjiity can be examined 
by considering Figure( 5.l6)3where the fluid parameters 
were allotred to vary. Since the added mass coefficient 
was computed for an equivalent annulus, its value is nrach 
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lower { = 2.125 as against^ = 4- in Figure 5»13). 

ihe dynamical behaviour is now quite different. Instabi- 
lities associated with first and second node are of ing 

type rather than oscillatory and occur at higher velocities 
= 1.83 and 3.66 respectively. Internal flow in sane 
and reversed direction ( ' Figures 5«17 and 5*18) reduc« 

these values to approximately 1.82 and 3-645 and 1 .81 and 
3*642 respectively. The third mode remains damped for the 
ra.nge of flow velocities considered. 

V.4 Cantilever Cylinder t 

The complex frequency diagram for the case of 
an isolated cantilever cylinder has been shovni in Figure 
(5*'I9)* The parameters correspond to those of Figure 6 of 
reference C55 3 • However, free end has been considered blunt 
and fluid parameters have been held constant. It can be 
seen that for loi; values of flow velocity the free motions 
of the cylinder are damped. As the flow velocity is gradua- 
lly increased the first mode locus crosses I’m(w)axis at 
”^2^' indicating the onset of flutter. The first mode 

loci bifurcate on Re ) axis, one branch. approaching the 
origin and the other receding ai-ray. The cylinder loses 
stability in the second mode by flutter, at xi 2 ^ 4.8. 

The third mode loci remain in the stable region for the 
range of velocity considered and inaginary parts of the 


frequencies of the third, node mono tonic ally increase vith 
flow velocity. 

The complex frecpaency diagram, shown in Figure 
(5*20) has "been obtained by allowing the fhaid parameters 
to vary. General character of the complex frequence?" diagram 
appears to be similar as that of Figure (5.19) except that 
the critical velocity for the onset of first mode flutter is 
now increased and subcritical damping is larger; the 
velocity for the second mode flutter to set in is, however, 
unaffected. 

The dynamical behaviour of the same cylinder \dien 
the direction of the flov7 is reversed, i.e. towards fixed 
end, can be examined from Figures (5.21) and (5.22) which 
correspond to figures (5.19) and (5.20) as far as other 
parameters are concerned. It can be seen that the loci of 
all mod.es leave the axis into unstable region, for 

small flow velocities. Tlius, flow induced damping becomes 
negative. It is iiiteresting to note that the first mode 
regains stability at U 2 - 0.4-1, _,but loses _ ' . 
stability by buckling at U2»^1.25. Second mode, similarly, 
regains stability at Ug^^^* 3-33. Third mode however remains 
unstable for the range of velocities considered. In case 
of cantilever cylinders subjected to steady flow tow’a.i-ds 
fixed end the flow appears to destabilize the cylinder for 
the flow velocities within the range considered here. 



For this case, vhen fluid, parameters have been alloved to 
vary, the results, shouri in Figure (5.22) were obtained. 

It is of interest that the onset of the first node flutter 
now disappears and first mode instability is of buckling 
type. The second mode loses stability by flutter at 
0.15 but regains it at u^«» 3.1 ; hov/ever, the third 
mode remains unstable for the velocities within the range 
considered for these results. 

The results obtained for cantilever cylinders are 
particularly striking; flutter Instabilities associated 
i"7ith cantilever cylinder are the result of frictional and 
inviscid hydrodynamic forces both of them being of non- 
conservative type in this case. Before examining the 
mechanism underlying these instabilities it will be con- 
venient to examine the effect of flows on the dj^namic 
behaviour of pipes . 

V.5 Plnned-Pinned Pipe ; 

Figure(5.23 )shows the complex frequencies of the 
lowest three modes of a pinned-pinned pipe with internal 
flovr velocity as parameter, for = 0.1, € = 'm-0, ^ = 1 

and ,o<s=K - V = TT = T = o . Pipe has been con- 

sidered to be immersed in an infinite fluid. Fluid 
parameters were computed as discussed in Chapter ilx, 
resulting in non— aero value for c,' the viscous drag 
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coefficient. It can be seen that small flow velocities 
damp the free motions and as the external fluid is sta- 
tionary, the damping, which is proportional to real part 
of the Complex frequency, is constant and roughly equal 
for all the three modes. The imaginary parts of the 
frequencies of all three modes decrease with the increase 
in velocity till those of first mode reach and bifurcate 
on the Re (w ) axis. One of the branches, then crosses 
the ajcis, at u^ , which indicates the onset 

of buckling. Further increase in velocity causes the two 
branches to receed away from the origin while the second 
mode- loci approach close to the Re ) axis. At 
u^ 6.08 two btanches of the first mode start approaching 
the origin. At slightly higher flow velocity second mode 
loci reach and bifurcate on Re ) axis. The pipe loses 
stability by buckling in the second mode at u^ssi 2 . 

With further increase in velocity two branches of the 
second mode receed away from origin till at u^ *5^ 6.39 
they coalesce with- the branches of the first mode loci 
and leave the Re ) axis where Re («h ) is positive and 
negative, indicating the onset of coupled mode flutter. 

In Figure 3(a) of reference C 55^ the coupled- 
mode flutter for the same case starts at u^«> 6.38, and 
has been reported there for the first tin©. FigureC5»24) 
shows the dynamical behaviour for this case when exter- 
nal flow velocity is 0.5) where the second mode loci 
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have been drawn on I>n(uj) axis for the sake of clarity, 
since they overlap those of the first njode. The flow 
induced damping for small flow velocities is increased 
due to external flow. General character of the frequency 
diagram is similar to the case 'vrhen external fluid was 
stationary, however, all the critical velocities are 
reduced. This observation is of practical interest since 
this case pertains to parallel flow heat exchanger. 

The complex frequency diagram of the lowest three 
modes for the pipe considered in Figure (5.23 ) (with the mass 
ratio ^ changed to 0,5) is shown in Figure (5.25). The 
presence of stationary fluid results in non-zero values 
for the parameters and c. The dynamic behaviour 

shown is different than the one shown in Figure 5.23. It 
is noted that once again the first mode frequency vanishes 
at u = 25 . Tt , However, u = 2 7X does not correspond to 
buckling in the second mode, but rather it is the point 
where the system loses stability by coupled mode flutter. 
It is interesting to note that similar results have been 
obtained in reference [ 48l for the case of an isolated 
pinned-pinned cylirder, for 0.4-8, = 0.25, 7 

=r = , Ihe role of the second 

mode in Figure (5-. 23 )is assumed in this case by the second 
branch of the first mode, without affecting the overall 
stability characteristics of the system in a radical way; 



this interchange in the roles of contiguous nodes has teen 
previously found to occur in the case of cantilever pipe 

C373 . 

The dynamic behaviour, when external velocity is 
0.5 in the same direction as that of the internal flow 
(refer to Figure 5.26) is not significantly different, tut 
critical velocities for buckling and coupled mode flutter 
are reduced, evidently due to frictional forces associated' 
with external flow. Eeversal in the direction of external 
flow (refer to Figure 5*27) does not cause significant 
change in the gener al behaviour excepting that the critical 
velocity for buckling is slightly reduced ( 3-09 as against 
3.11 in Figure 5*26). 

V,6 Glamned-Clamped Pipe 

Finally the case of a pipe with clamped-claniped 
ends may be examined. The results are shown in Figure 
( 5 . 28 ), corresponding to parameters T - TT ~oL^ 

P = 0,5 and (2 = ^0. It is seen that the general character 
of the results is similar to those obtained for a pi.ined- 
pinned pipe, shown in Figure (5 •25.). The velocities at 
which system loses stability • ■ by buckling and common 
mode flutter are higher due to the nature of end conditions 



V.7 Comparison Between Cylinders and Pj-pes 


The similarities and differences in the results 

for pipes and cylinders with different end conditions 

can now he examined together in general terms. General 

behaviour of the complex freauencies with increasing flow 
-pinned 

for pinned/pipes and cylinders is similar except that for 
pipes second mode flutter does not set in, at least for 
the range of parameters of this study. The difference 
between systems of high and low mass ratio is worth 
considering. For sraall values of ^ the system buckles 
in the first mode both for pipes and cylinders; but the 
second mode loses stability by flutter in case of cylinder 
and by buckling in the case of pipe. The loss of stability 
by coupled-mode flutter is common to pipe and cylinder. 

For larger values of ^ , second mode instability is 
notably absent in the caae of a pipe while in the -ase of 
a cylinder, common mode flutter is absent. The coupled 
mode flutter was also found to set in for thin pipes 
studied by Paidoussls [431 , Weaver and Unny [49] and 
weaver and Mfklatun [ 50 1 . The results for clamped- 
clamped pipes are similar to those for pinned-pinned 
cylinders. For supported systems flow induced damping 
at small flow velocities is caused only by external flow. 
The dynamical behaviour discussed above also holds for 
the cylinder or pipe in a cluster but higher added mass 



The - character 


coefj. icient dsetahilizes the systeia |_4-83 
of the instahilities highly depends on the added mass, 
as shown "by the results with low and high added mass 
coefficients. Thus the geometry and spacing between the 
cylinders are parameters of equal importance in clusters. 

The effect of flow on the free motions of canti- 
lever cylinder results in oscillator;/ instabilities of 
all the three modes studied. This can be attributed to 
the irictional and Inviscid hydrodynamiic forces, as the 
system is iniierently non conservative. The effect of 
two flows in all the cases considered seems to destabilize 
the system, the effectiveness of destabilization being 
dependent on the relative direction of the flows, for sys- 
tems with symmetric supports , 

V.8 The Mechanism of Instability 

Having seen the nunBrical results obtained for 
cylinders and pipes subjected to internal and external 
flows the mechanism causing various instabilities may now 
be considered from the analytical point of view. 

V . 8 . 1 Buckling Instabilities 

The mechanism underlying buckling instabilities 
may be shown by considering the static equilibrium of the 
cylinder with both ends supported, which is assuraed to 



have momentarily frozen taking an arbitrary bovred shape 
j (x). Eliminating the ticje-dependent terms in equation 
(2.20) and considering excitation f (x, t) to be absent, 
one obtains 


oi. 4 _q 


y 




-b '^-2.0 


4^ 

cix.^ 


- O 


(5.1) 
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- 5^ w,g^- A r If q-P, )-b?ybyA,/i>*0] 

+-i < '-i) -^1 

Cn- 3) ? itObkCb 


For sufficiently high fluid velocity the forces 
represented by the terms and ^ 

may overcome the flexural and tensila restoring forces 
P T ^ ^ and -5T ^r-^, respectively, resulting 

Irp- 

in a monotonic increase in the amplitude of y (x). It 


niay further be noted that the centrifugal forces due to 
inner and outer flow add irrespective of their direc- 
tions, resulting in lower critical velocity for buckling 
when both the flows are present. In case of pipes 
carrying fluid, when gravity is ignored, equation (2.21 ) 
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sirnplifies to (with time dependent terns eliminated 


), 


where 


h ^0- 


cl\ 




-O 


(5.2) 


u"' - & L(i-2."i») n + r] 

The above equation admits solutions of the tyre 

> lY-?* 

i A-e 

3-' 

mth = V-2_- o 5, Y3 - “'O'? T4. . Hence the solution 

can he rewritten in the form 

y\i % ) = A , + A2.e + 6,^°^ 

Substituting the boundary conditions (equations (2,22), 
for simply supported pipes) for non- trivial solution one 
gets a vanishing determinant which yields 

Sn\ \> =. o 


Since “O' = 0 is a trivial solution, one obtains 


(5.-) 


where j = 1 5 2, 


This, of course 5 is the familiar 


Euler’s result if pipe is considered as an equivalent 

column subjected to a compressive load j^AUbj 

thus in case of u =J\.= 0, one has , the 


compressive load at the end of the pipe for the onset 
of buckling. In case of clamped-clamped pipes critical 
values of are found corresponding to the zeros of 

Z ( I- (L05\^l - Vi'" ^ = “ 

Ehysically, one may regard 'iT ss an effective 
compressive load. Alternatively, ano. mere appropriate!^' 
for the case of a flo'vring fluid, H may rtjgar.-ei^ 


(5.5) 



as a generalized centrifugal force . When ttiis force 


overcomes the flexural restoring force, the pipe huchles. 


V.8.2 Flutter Instabilitie s 


For oscillatory instahilities, the condition of 
neutral stability is one of the dynamic equilibrium where, 
in the course of one cycle of oscillatioiij the eneXvgy 
transfer from fluid to cylinder and vice versa exactly 
balances. When the former exceeds the latter (since the 
fluid stream may be regarded as a source of infiriite 
energy) the amplitude increases without limit. In the 
opposite case, the oscillations are damped. 


Considering equation (2.19 a)? one finds that 
the rate of work done on the cylinder in the course of 

free periodic motions may generally be written as 

U _ U _ ^ . . 2 . 


d W 


d-t 

u 


- ( El 3 \ ^Adu*- 


! B :>c / 




S '5( 


■b 'JbJb.j.I ') 

and over one period of oscillation t-| , the work done 
, with the cylinder returnJ.ng to its original 
position is given by 


( 5 . 6 ) 


f,A Ub^b.C y 


AW = - \ ^ 



•^pplyi^g "the boundp.ry conditions, one obtains for 


cantilever cylinders 
■fc ! 


r ^ ^ i i ^ i5> Ji ) 


t-t L 

rV 


1 J C Ut>2_St:5^ ^ )) ) dx^it 


O O 


( 5 . 8 ) 


( 5 . 9 ) 


and for cylinders with both e^ds supported 

aw= ”4 ^ )dxdt 

o o 

If AW Co , motion will be damped | wiiile if C.W > O' 
the oscillations will grow, i.e. the s^’^stem will be unsta- 
ble by flutter [ 29 ] . From equation (5.2) it is evident 
that in the absence of hydrodynamic drag effects, in case 
of supported cylinders, AW = 0. Thus the hydrodynamic 
forces are of purely conservative type. It is, therefore, 
evident that the existence of fliitter for supported 
cylinders depends on the presence of frictional forces as 
the system is inherently conse]r.’'ative . This explains 
the onset of' instability bj'- ordinary flutter in case of 
supported’ cylinders and pipes when the external flow is 
present, on the other hand, cantilever cylinders iray 
become unstable entirely by the action of inviscid 
hydrodynamic forces. It is not surprising tnat a small 
correction in the viscous terms of the equation of motion 
should have a profound effect on the dynamical behaviour 
of cylinders where the system is iiAierentlj^ conservative. 
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The onset of amplified oscillations corresponds 
to a cnange in the character of the mode at a specified 
flow “velocity as determined "by the conseri/ative forces 
in the s^^'ctenij frictional effects, however small, cause 
the amplification of the oscillation, whereas before the 
change oc cured they caused damping. Tiriis change in the 
character of the mode is essentially like the change from 
class B to class A behaviour discussed by Benjamin L15 1* 

Inspecting equation (5*9) more closely, one 
observes that the condition of neutral stabilit^r (CiW =0) 
requires that 

o 

For the right hand side of equation (J.10) to be nonzero, 
different parts of the cylinder must vibrate in quadrature, 
and this is an effect that can only be produced by the 
action of the flow. If the solution is represented as the 
sum of a standing and a travelling wave, it can be seen 
that the right hand side of equation (5.10) depends only 
on the amplitude of the travelling wave component 5 further- 
more , the phase velocity of this component must be in the 
direction of the flow for the right hand sine of equation 
(5.10) to be positive. 

In the case cf cantilever cylinders, from equation 
(5.8), it can be observed that the energy transfer between 




( 5 . 10 ) 



the fliiid stream and the cylinder may he caused loth by 
inviscid and viscous hydrodynamic forces. EGnja.!i;in ^ 11^ 
has considered the mechanism of ener^ry transfer in the 
related problem of fluid flouung in an articulated pipe 
system. He found that over one complete period of 


oscillation t^ , the energy gained by pipes when the doT,.ni 
stream end is free, is given by 

AW-- MiU (5.11) 

o 

where is the mass per unit length of the contained 

fluid and T and R are the tangential and position 
vectors at the end of the last pipe. Similarly between 
equation (5.1 1) and the first two terms of equation (5*6) 
is striking. In the case of cantilever cylinders, the 
oscillatory instabilities can occur independently of the 
friction forces. T/iliien = 0, the equation (5*8) gives 

A 'V = - f V, A U&fci vt t )1 

D • , 1, 

o 

It is evident that when \ is small, the first 
term in each integrand predominates ovei' tne second and 
vibrations are damped. For sufficiently high fluiu 
velocity^ however, if y and y^ arc sufficiently out of 
phase to give a. negative average value to y y , it is 


( 5 . 12 ) 



possible for to become positive, in ■which case 

oscillations are amplified. Thus for amplified vibration, 
for the greater part of a cycle the free end of the cylin- 
der must slope baclowards to its motion, i.e, a dragging 
motion must occur, 1 as given by equation (5.12) is the 
work done by nonconservative hydrodynamic forces. Equation 
(5*12) also shows that for the oscillatory instability the 
directions of the two flows are important in case of 
cantilever C 3 ?'linders. This is particularly clear from the 
numerical results. 



The coupled-mode flutter at sufficiently high flow 
velocity was shown to occur in the case of thin pipes 
conveying fluid. The dynamics of this was studied using 
shell theory to describe the motion of the pipe ^ 4-3, ^6, 

4-9, 50~^ . The coupled mode, flutter was obtained for thicker 
pipes only recently \ 55 I • The existence of the coupled- 
mode oscillatory instabilities is connected to the presence 
of Coriolis forces which are operative only when mass 
ratio [3 ^ 2. He confirmed this by calcula'cing tne fre- 
quencies for a simply supported pipe with {- =0. Due 

to the effect of Coriolis forces tiie '’critical fj.ow 
velocities obtained by Euler’s method of equilibrium do 
not always correspond to thresholds of, instabilities, as 



96 


indicated by the results obtained using dynamic method. 

The systems with supported ends considered in the present 
study are subject to inactive Coriolis forces and there- 
fore, are gyroscopic conserva.tive systems. Accordingly, 
the static method gives only the first critical velocity 
which corresponds to buckling in the first mode; beyond 
this point systems cease to be positive definite, and 
only the dynamic method is capable of giving all higher 
critical velocities and their correct physical interpre- 
tation. From equation (5*9) it is seen that in the absence 
of frictional forces, /\y=:o conservative systems 

considered here. On closer examination, this argument 
implies that the critical point for flutter cannot be a 
point of neutral stability. The critical condition for 
coupled-i;iode flutter is not a point of neutral stability 
as it corresponds to coincidence of two frequencies and 
hence to growing oscillations of the form y ■■ , tb = 

^ ( xt C cv + b t ) i , , Therefore one may arrive at 

a conclusion that, apart from buckling in the absence of 
frictional forces the gyroscopic conservative sj^s terns can 
only be subject to coupled mode flutter. The presence of 
ordinary flutter for these systems then can only be 

connected to the hydrodynamic drag forces. Recently 

Shieh ! 6^] has drawn similar conclusions for gyroscopic 
conservative systems in general. 



Ihe existence of Coriolis force in the systeE 
gives rise to phase difference and mode coupling. Coriolis 
force terms introduce asymmetrjr in damping matrix C anc?. 
inertia me,trix K of equation (2.28) and make the problem 
non-self-adjolnt . As shovm 83’' Chen (,4-01 , Coriolis forces 
do not dissipate or supply any energy?". Coriolis forces 
contribute to de.mping only when displacement a.t one of tlie 
ends ot the cylinder is nonzero. For simpli' supported pipe 

or cylinder, from equation (4-,l4-), 

1 

G • = 0 for n, 3 = even or n, j = Odd 

1J.J 

This relationship implies that there is no coupling between 
even-even modes or odd-odd modes. Physically this means 
that the motion in odd modes v/ill excite the motion in 
even modes only, and vice versa. For the fundamental mode 
of pipes or cylinders with supported ends, this can be 
explained as follows: the displacement being symmetric 
about the centre of the span, the angular velocities of 
left and right half span of the beam are in opposite direc- 
tions. The fluid motion being unidirectional, gives rise 
to the antisymmetric Coriolis force which excites only 
antisymmetric modes. The motion of these modes in turn 
produces a symmetric Coriolis force which induces sj^'inme— 
trie motion. This is the reason wh3’' Coriolis force causes 
mode coupling and phase difference. This further explains 
to some extent, the mechanism giving rise to coupled mode 
flutter at higher flow velocities. 
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CmPTER VI 


COIJCLTJSIONS 

VI. 1 Sumr.ary rnd Results 

In this chapter the results obtained in Chapter 7 
have been suinir;arized and the 3 " are reexaiiiined in the con- 
text of theory developed in Chapter II and objectives set 
out in article 1,4-. 



As a first step in understanding the mechaniSE of 
parallel flow induced vibration, mat!'.ematical model 
governing small transverse vibrations of c^ainders 
conveying fluid and subjected to external steady flow 
has been set up taking into account gravity, pressuri- 
zation, internal dissipa.tion and direction of the 
flows. The equation of motion as derived in the present 
study corrects an error in the representation of fric- 
tional forces, considered in previous formulations 
r^-s, 6 i 3. Galerkin method has been used to obta.in the 

■ . I 

complex frequencies of the cylinders' vr^h different 
end conditions. The orthonormal eigen functions of a 
simple beam with identical boundary conditions in the 
absence of fluid have been used as coordinate func- 
tions for Galerkin’s expansion. 
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2. In the case of pinned-pim^d (siEiply. supported) 
cylinders and pipes small flow velocities tend to 
damp the free motion, as the flow velocity is gradually 
increased, the cylinders and pipes are subject to 
buckling. As velocity is increased further, cylinders 
lose stability by flutter while in the case of nipes 
the second instability’" is by buckling in the second 
r.ode if the mass ratio is small and by coupled mode 
flutter if mass ratio is high. With still higher flow 
velocity, coupled-mode flutter has been found to occur 
in the case of cylinders as well as pipes, especially 
when mass ratio is hi^. 

3* When the cylinder or pipe is a part of a cluster of 
identical cylinders or pipes the close spacing results 
in increased virtual ciass, which has been found to act 
as destabilizing agent. Further, for clusters, in 
addition to ordinary flutter the coupled-mode flutter 
has been found to occur t^^ice with buckling instability 
interposed in betwen. When the "virtual mass coefficient 
is coKiparatively small, on the other hand, only buckling 
instabilities in the first and second mcxie have been 
obtained for sufficiently high flow velocities. 

Similarity in the character of "these ins "tabili ties and 
those in the case of supported pipes with low mass 
ratio maybe noted. 



t 

4. The presence of second flow (internal flow in case of 
cjainders and external flow in case of pipes) has teen 
found to loi/er the critical velocities in general. 

The severity of the destabilization effect bein^, depen- 
dent on the directions of the two flows. The effect on 
critical velocities for buckling is less compared to tha 
on the velocities for the onset of flutter. 

5 . In case of cantilever cylinders subjected to only ex- 
ternal flow toxirards the free end, small flovi velocities 
have been found to damp the free raotion. However, the 
flutter involving first mode soon sets in, -followed 

by one involving second mode. When the direction of 
the external flow is towards fixed end, the flow induces 
negative damping for small velocities. As velocity is 
increased further the cylinder regains stability 
Cu2«5s 0.41 ) and looses it again (u 2 ^ 1.25) by diver- 
gence involving first mode. Similarly, second mode 
regains stability, at still higher velocity (U 2 ®* 3*33); 
however third mode frequencies remain in the unstable 
region for the range of velocities under consideration. 

6. 'The divergence type of instability has been shown to 
be due to inviscid centrifugal force of the fluid, 
which overcomes the elastic and tensile restoring 
forces at sufficiently high flow velocities. This 
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explains destabilization effect of added mass in ce.se 
of clusters and that due to the presence of both in- 
ternal and external flovs.in the case of pi pcS £.nci 
cylinders. The direction of. the internal flow dees 
not affect the critical velocity for buckling. 

7 . Instability by ordinar]/- flutter is caused by \'iscous 
forces due to external flow alone in the case of suppor- 
ted systems and by inviscid hydrodynamic forces in 
addition to viscous forces in the case of cantilever 
cylinders and pipes. For oscillatory instabilities tlie 
direction of the external flow is important when inter- 
nal flow is also present. The correction in the visccais 
terms of equation of notion has resulted in ordinary 
flutter instability in case of cyliMers and pipes mth 
external flow. In case of cantilever cylinder, which 
constitutes an inherently nonconservative system, ttiC 
oscillatory instabilities are due to hydrodynamic 
forces. 

8. The coupled mode flutter has been explained to be the 
result of inactive Coriolis force which introduces mode 
coupling and phase difference. The static method to 
obtain critical velocities for buckling in higher modes 
fails due to the presence of gj^roscopic Coriolis force, 
the system no longer remaining positive definite. 
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In the parameter space, the critical point for the 
flutter cannot be a point .of neutral stability. In 
the absence of frictional forces supported pipes and 
cylinders can onl.y be subject to coupled-mode flutter 
since the systems are conservative. 

VI. 2 Future Work ; 

1 . The results concluded so far constitute the first 
step towards the solution of response to the flow 
induced excitation which may consist of either fluc- 
tuations in the wall pressure due to turbulence or 
parametric excitation due to piasations in the flow 
rate as discussed in Chapters I aiil Illi The condi- 
tions , under which the bouixied solution for the motion 
cannot be obtained, have only been exaroined in the 
present work. The next step is to seek the solution 
of the response problem. If the excitation is para- 
metric, the stability problem has to be solved for 
unsteady flow velocity. In particular, when velocity 
has a SBKtll haraonic perturbation superimi)osed on a 
steady mean value, the primary and secondary instabi- 
lity regions need to be investigated in the context 
of corrected viscous interaction and presence of 
internal ard external flows. 
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2. It I'ay te noted that though the methematicai Eodel 
formulated in the present work takes an account of 
gravity, internal dissination, pressurization and 
spacing tetween the cylinders in a cluster, detailed 
paranetric study has not been taken up. Further \'ork 
on this topic should take up the study of the effect 
of these parameters, in addition to the structural 
and fluid coupling between the cylinders, on the dyna- 
mic behaviour of the cluster as a whole, 

3- The present study, as mentioned in Chapter I, had to 
be diverted and restricted to the investigation' of 
stability anal3/'sis only due to the non-availability 
of reliable analytical or experimental data. This is 
due to difficulties in the measurement of snail scale 
pressure fluctuations in complicated flow geometries 
of practical interest. Analytically, the pressure 
fluctuations can be obtained as mentioned earlier but 
the method is uneconomical; furthermore it would be 
more convenient to obtain the response of the cylinder 
directly using simulation technique rather than pre- 
ssure fluctuations. It may be a-ppreciated that if 
the excite, tion contains de to mini 3 tic unsteady forces 
or force field in addition to random but stationary 
forcing function, as in turbulence induced excitation 
due to pulsatile flow of a fluid, no satisfactory* 



rret::cd c:.:ists for obtaining tho response. Consid:ri:i£: a 
these difficulties inanediato need of approxiffiate liXti-ods 
for obtaining the response exists. This is rise one of 
the reasons for largo volurc of research acti-'/itp on sta 
bility probleffi conpered to the one on response problen. 
An alternative offered by experimental measure- ■•onts of 
the response itself, though useful for particular system 
is not sufficiently attractive due to difficulties in 
extrapolating the results. Under the circumstances the 
measurement of pressure fluctuations and simlation are 
the only feasible methods. In theory, one can obtain 
the response coupled vith flow dynamics, even when flow 
velocity contains unsteady components. This \iray, even 
though the solution is obtained by computationally un- 
economical methods, one has a satisfaction of obtaining 
a meaningful solution. 
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